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EXECUTIVE  SUMMARY 


MAGLEV  trains,  electromagnetic  railguns  and  electric  rollercoasters  are  all  examples  of 
linear  motors,  a  technology  that,  through  years  of  development,  is  finally  coming  of  age. 
A  cousin  of  the  linear  motor,  the  linear  generator  has  historically  not  received  the  same 
level  of  focus  and  development.  The  linear  generator’s  level  of  attention  is  growing, 
however,  and  new  and  emerging  technologies  that  incorporate  linear  generators  are 
becoming  less  niche  and  more  mainstream. 

Wave  energy  converters,  including  wave  surge  converters,  surface  and  subsurface 
point  absorbers,  and  surface  attenuators,  all  stand  to  benefit  from  the  inclusion  of  linear 
generators  in  their  existing  designs.  Some  companies  have  begun  to  develop  and  test 
wave  energy  converters  that  include  linear  generators,  but  the  vast  majority  of  wave 
energy  production  is  still  done  through  hydraulic  turbine  technology.  The  Department  of 
Defense  and  the  United  States  Navy  stand  to  benefit  from  the  development  of  this 
technology  for  use  in  small-scale  wave  energy  production  to  power  remote  sensors  and 
recharge  unmanned  vehicles  without  the  need  of  a  support  vessel. 

Free-piston  devices  are  another  example  of  a  burgeoning  use  of  linear  generators. 
These  devices  can  be  either  Stirling  engines  or  more  conventional  internal  combustion 
engines.  In  both  cases  the  reciprocating  piston  is  driving  a  linear  generator,  resulting  in 
direct  electrical  power  production  and  eliminating  several  energy  conversion  processes 
along  the  way.  Free-piston  engine  designs  also  naturally  lend  themselves  to  the  use  of 
different  fuels  from  gases  like  hydrogen  and  natural  gas  to  combustible  fluids  like 
kerosene  and  gasoline.  These  devices  are  ideal  as  electric  vehicle  range  extenders,  and 
they  are  already  in  development  by  major  brands  like  Toyota.  This  flexible,  modular 
design  also  lends  itself  to  distributed  backup  electrical  generation  on  the  Navy’s  electric 
ships  of  tomorrow. 

Solid  state  power  electronics  is  an  enabling  technology  for  linear  generator 
development.  Finear  generators  by  their  nature  come  to  a  stop  twice  per  cycle,  when  the 
translator  reaches  either  end  of  its  travel.  At  this  point,  the  generator  is  producing  no 
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power.  This  start/stop  action  causes  a  less-than-ideal  power  quality  on  the  output  bus. 
This  can  be  seen  in  the  simulation  results  as  well.  It  is  thanks  to  advances  in  power 
electronics  that  the  output  of  linear  generators  can  be  transformed  into  usable  power. 

Mathematically,  a  linear  generator  is  very  closely  related  to  a  standard  rotary 
generator.  For  this  reason  the  generator  model  is  based  on  a  rotary  machine  described  by 
Chapman  et  al.  in  [1].  This  machine  is  a  three-phase,  non-salient,  non-sinusoidal 
permanent  magnet  synchronous  motor.  The  variables  in  the  model  used  for  this  thesis  still 
reference  the  angle  and  angular  velocity  of  the  rotor  despite  the  linear  generator  having  a 
translator  that  is  described  in  terms  of  speed,  position  and  force.  External  to  the  model, 
the  torque  of  the  rotor  is  converted  to  the  force  of  the  translator;  similarly  the  angular 
velocity  is  converted  to  speed  and  the  rotor  angle  to  translator  position.  The  model  also 
employs  a  rotor  reference  frame  using  a  qdO  transformation,  a  method  closely  related  to  a 
Park’s  transformation. 

The  advantages  of  using  a  qdO  transform  do  not  apply  to  this  machine  because  it 
has  a  non-sinusoidal  back  electromotive  force  (EMF),  as  seen  in  Figure  1.  The 
transformation  is  included,  however,  for  added  utility  in  future  work.  If  the  values  of  the 
harmonics  in  the  flux  distribution  are  set  to  zero  then  the  model  becomes  sinusoidal  and 
the  advantages  offered  by  the  qdO  transform  of  simplified  control  and  constant  state 
variables  in  steady  state  conditions  can  be  utilized. 

For  the  simulation,  each  phase  of  the  generator  model  feeds  a  full  bridge  diode 
rectifier.  The  outputs  of  the  three  diode  rectifiers  feed  a  common  dc  bus  and  a  150Q  load 
resistor.  Only  the  bus  voltage  is  shown  because  the  load  is  purely  resistive. 
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Time  (arbitrary) 

Figure  1.  Non-sinusoidal  back  EMF. 

Several  different  mechanical  torques  are  used  as  inputs  to  the  simulation.  In  the 
case  where  a  wave  energy  converter  is  being  emulated,  the  simulation  uses  a  IHz  low- 
amplitude  sine -wave  mechanical  torque  as  the  input  function.  This  input  torque  is 
presented  as  the  equivalent  mechanical  force  on  the  translator,  as  seen  in  Figure  2.  The 
linear  generator  does  produce  a  useful  amount  of  output  power  in  this  scenario,  but  the 
power  is  of  poor  quality,  as  seen  in  Figure  3.  The  output  voltage  has  a  ripple  that  is 
greater  than  half  of  the  bus  voltage,  which  can  be  seen  by  comparing  Figures  2  and  3  as  a 
function  of  the  absolute  value  of  the  input  force. 


To  use  this  type  of  power  as  an  ac  source  or  a  regulated  dc  source  more  work 
must  be  done,  either  with  filter  components  or  with  active  power  electronics.  If  the  power 
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into  the  machine  (torque  times  angular  velocity)  is  calculated  and  compared  to  the  power 
out  of  the  machine  (voltage  times  current)  this  scenario  sees  the  model  generator 
operating  at  54.1%  efficiency. 


In  a  different  setup,  the  simulation  uses  an  impulse  that  fires  additively  in  the 
direction  of  travel  of  the  translator  each  time  the  translator  passes  the  midpoint  of  the 
generator.  The  simulated  response  of  the  model  generator  can  be  seen  in  Figure  4,  where 
the  parameters  are  scaled  in  order  to  best  be  displayed  together  in  an  intuitive  manner. 


Figure  4.  Input  and  spring  forces  acting  on  the  translator  and  its  resultant  motion. 
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The  Fspring  in  Figure  4  represents  the  foree  of  a  eentering  spring  that  is  attached  to 
the  translator  to  aid  in  the  stability  of  the  system  and  in  this  case,  as  a  return  force  for  the 
translator.  This  setup  draws  parallels  to  the  free-piston  engine  arrangement  already 
discussed;  most  free-piston  devices  employ  an  air  cushion  spring  to  act  as  a  return  force. 

The  frequency  of  the  input  force  in  this  scenario  is  controlled  by  behavior  of  the 
machine.  In  this  simulation  the  machine  stabilized  at  a  frequency  of  70.4Hz  or  about 
4225RPMS  if  equated  to  the  terms  of  an  internal  combustion  engine  (although  there  is  no 
rotation  here).  This  higher  frequency  results  in  a  marked  improvement  in  power  quality 
on  the  dc  bus.  This  can  be  seen  in  Figure  5,  where  the  voltage  ripple  is  reduced  to  2.7% 
of  the  total  dc  bus  voltage.  The  machine  also  exhibits  an  improved  efficiency  of  71.2%. 
This  improved  efficiency  is  somewhat  expected  since  it  is  in  keeping  with  the  trends  laid 
out  by  Eric  Laithwaite’s  goodness  factor  formula,  which  states  that  machine  efficiency  is 
proportional  to  its  speed  among  other  factors. 


Figure  5.  Output  voltage  on  the  dc  bus  from  a  floating  frequency  impulse  input 

torque. 

Through  various  simulations  the  Simulink  model  proved  to  be  a  useful  and 
flexible  tool  in  analyzing  the  performance  and  operating  characteristics  of  a  linear 
generator. 
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There  is  more  work  needed,  to  include  some  of  the  non-ideal  characteristics  of  a 
generator  such  as  demagnetization  and  core  saturation.  The  model  could  also  benefit 
from  the  inclusion  of  cogging  torque,  especially  in  the  low  speed  situations  where  its 
effects  would  be  the  most  observable. 
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I.  INTRODUCTION:  A  BRIEF  OVERVIEW  OF  ELECTRICAL 
MACHINES,  THEIR  ORIGINS  AND  USES 

A,  BACKGROUND 

In  1820,  Hans  Christian  0rsted  discovered  he  eould  make  a  eompass  needle 
deflect  when  he  passed  eurrent  through  a  nearby  wire,  thus  uncovering  a  relationship 
between  eleetrieity  and  magnetism.  Just  eight  years  later,  in  1828,  Anyos  Jedlik  had 
already  ereated  the  flrst  eommutated  rotating  eleetromeehanioal  maehine.  In  1831, 
Miehael  Faraday  diseovered  that  a  eurrent  eould  be  generated  by  varying  a  magnetie 
fleld,  diseovering  induetion  and  paving  the  road  for  Hippolyte  Pixii’s  invention  of  the 
eleetrie  generator  a  year  later.  Shortly  thereafter,  in  1837,  the  flrst  United  States  patent 
was  awarded  for  an  eleetrie  motor  to  Thomas  Davenport,  it  was  U.S.  Patent  Number  132 
[1],  [2],  [3]. 

The  first  U.S.  patent  for  a  linear  maehine  would  not  oeeur  for  another  68  years;  it 
would  be  awarded  to  Alfred  Zehden  under  the  deseription  of  “Eleetrie  Traction 
Apparatus”  in  1905  [4].  Zehden’s  invention  was  the  forerunner  to  the  modern  day  maglev 
system.  One  of  several  eonfigurations  deseribed  by  Zehden,  in  his  patent,  is  shown  in 
Figure  1. 


Figure  1 .  Finear  motor  propulsion  system  as  shown  in  U.S.  Patent  No.  132, 

from  [4]. 
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Progress  in  linear  motors  has  been  greatly  fueled  by  interest  in  maglev  technology 
for  transportation  and  by  military  interest  in  devices  such  as  the  electromagnetic  railgun 
and  the  electromagnetic  aircraft  launch  system.  Recently,  there  has  been  a  growing  and 
significant  application  of  linear  motors  in  industry  for  use  in  robotics.  Linear  generators, 
on  the  other  hand,  have  had  their  progress  chiefly  fueled  by  interest  in  new  forms  of 
energy  production  which  will  be  discussed  in  detail  in  Chapter  II. 

An  enabling  technology  for  linear  generators  is  the  advancement  of  permanent 
magnets.  While  some  of  the  earliest  electric  machines  employed  the  use  of  permanent 
magnets,  those  machines’  capacity  and  usefulness  were  severely  stunted  by  the  lack  of 
strength  that  early  magnets  suffered  from.  Developed  in  1982,  the  neodymium  magnet 
can  achieve  a  remanence  (Br)  of  more  than  1.6  teslas  (T),  with  1.3T  being  representative 
of  their  standard  strength  [5].  It  is  true  that  neodymium  magnets  are  stronger  than  other 
forms  of  permanent  magnets,  but  also  that  they  are  more  expensive  by  weight  than  most 
forms  of  permanent  magnet.  To  understand  why  neodymium  magnets  are  utilized  to  such 
a  large  extent  one  must  compare  cost  per  maximum  energy  product  (BHmax),  in  units  of 
megagauss-oersteds.  BH^ax  is  also  referred  to  as  magnet  performance  and  is  often  used  as 
a  rating  system  for  magnets.  A  comparison  of  different  magnet  materials  and  their 
relative  costs,  by  magnet  manufacturer  Integrated  Magnets,  can  be  seen  in  Table  1  where 
NdFeB  is  the  neodymium  material. 


Table  1.  Relative  costs  of  different  permanent  magnet  materials,  after  [5]. 


Material 

BHmax 

(MGOe) 

Relative  Cost 
{$/lb) 

Relative  Cost 
{S/BH„,3x) 

NdFeB 

40 

35 

1.7 

SmCo 

26 

60 

4.9 

Alnico 

5 

25 

9.5 

Ceramic 

3 

2 

0.9 

Flexible 

1 

1 

1 

Another  important  factor  to  consider,  when  discussing  the  economics  of 

neodymium  magnets,  is  the  added  benefit  of  miniaturization  of  the  devices  the  magnets 

are  used  in.  Using  neodymium  magnets  allows  electric  machines  to  have  higher  force 
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densities,  allowing  for  smaller  and  eheaper  maehines  or  more  powerful  maehines  for  a 
given  size  depending  on  the  constraints.  The  effect  that  magnet  strength  has  on  the 
performance  of  an  electric  machine  will  be  discussed  in  detail  in  Chapter  III,  but  it  is 
evident  that  the  development  of  neodymium  magnets  has  had  a  profound  effect  on  the 
performance  of  permanent  magnet  electrical  machines  and  has  enabled  new  and 
compelling  applications  for  them. 

B.  POTENTIAL  INTEREST  TO  THE  DEPARTMENT  OF  DEFENSE  AND 

U.S.  GOVERNMENT 

The  Department  of  Defense  (DOD)  has  been  spending  a  lot  of  money  developing 
linear  motors  in  recent  years,  agreeing  to  pay  $676.2  million  for  the  development  of  the 
electromagnetic  aircraft  launch  system  (EMALS)  for  use  on  Ford  class  aircraft  carriers 
and  another  $240  million  on  phase  one  of  the  U.S.  Navy  electromagnetic  railgun  system 
[6],  [7].  Linear  generators  do  not  share  this  level  of  funding  or  attention  but,  with 
growing  interest  in  the  government  in  green  energy  sources  and  distributed  electricity 
generation,  some  might  argue  they  should. 

1,  The  Electric  Ship  and  the  Future  of  Electricity  Generation 

Emerging  technology,  such  as  the  pairing  of  a  free-piston  crankless  internal 
combustion  engine  with  a  linear  generator,  allows  for  a  more  direct  production  of 
electricity  than  traditional  means.  In  this  configuration  the  piston  becomes  the  translator 
of  a  linear  generator  with  the  cylinder  acting  as  the  stator.  One  advantage  of  this  setup  is 
that  it  reduces  mechanical  losses  caused  by  gearing.  It  also  has  the  potential  to  eliminate 
drivetrain  parts  such  as  transmissions,  reduction  gears  and  separate  electric  generators. 
Eliminating  portions  of  the  drivetrain  will  reduce  weight  and  cost  and  will  improve 
reliability. 

This  setup  is  also  very  suitable  to  a  modular  design,  allowing  for  flexible 
electricity  generation  capacity  to  meet  the  demands  of  separate  mission  packages.  The 
nature  of  a  standalone  modular  design  would  also  allow  for  distributed  generation  across 
a  ship  leading  to  better  survivability.  While  traditional  engines  used  in  direct  drives,  such 
as  marine  diesel  or  gas  turbine,  are  less  efficient  at  lower  power  settings  where  they 
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operate  for  the  majority  of  the  time,  the  modular  free-piston  engine  can  always  runs  at 
peak  efficiency  by  changing  the  number  of  pistons  operating  instead  of  the  load  of  each 
piston  in  order  to  meet  demand. 

The  simplicity  of  this  system  has  led  to  NASA  investigating  a  related  free-piston 
Stirling  engine  coupled  to  a  linear  generator  for  use  in  space.  This  is  discussed  more  in 
Chapter  II  but  is  mentioned  here  to  emphasize  the  simplicity  and  reliability  that  such  a 
system  promises. 

2,  Unmanned  Sensor  Energy  Needs 

Whether  they  are  stationary  sensor  arrays  or  autonomous  underwater  vehicles 
(AUVs),  the  U.S.  Navy  has  an  interest  in  sensing  more  areas  of  the  ocean  utilizing  less 
manpower  and  resources.  The  current  fleet  of  AUVs  operated  by  the  Navy  is  primarily 
powered  by  batteries  which  require  the  vehicles  to  be  recovered  every  few  days  or  even 
every  few  hours  so  they  can  be  recharged  [8].  There  has  been  research  on  AUVs  using 
inductive  charging  to  recharge  their  batteries  underwater  utilizing  shore  or  surface  based 
power  transmitted  via  undersea  cables.  Undersea  power  sources  are  non-existent  in 
remote  areas  of  the  ocean,  however,  so  in  order  to  reduce  the  resource  requirements  and 
increase  the  capabilities  of  the  AUVs  an  independent  power  source  must  be  developed. 

The  focus  of  the  Navy  is  increasingly  toward  the  littorals  where  shallow  water 
depth  presents  an  obstacle  to  a  fleet  that  is  designed  to  be  effective  in  deep  water.  AUVs 
are  a  big  portion  of  the  emerging  solution  to  this  dilemma  since  they  operate  effectively 
in  very  shallow  water  depths,  but  their  support  vessels  cannot  always  do  the  same.  This 
same  shallow  water  is  ideal  for  many  designs  of  wave  energy  converters.  While  the 
majority  of  the  development  in  wave  energy  converters  is  on  a  commercial  scale,  it  is 
feasible  to  scale  them  down  to  cheaper,  simpler  versions  of  existing  designs.  For 
example,  DARPA  worked  with  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI) 
to  develop  a  power  generation  buoy  in  the  Monterey  Bay  that  developed  more  than  400 
watts  in  normal  conditions  [9].  This  buoy  utilized  hydraulic  pressure  to  generate  power, 
but  it  is  a  perfect  example  of  where  a  linear  generator  could  have  been  substituted  to 
simplify  the  design. 
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The  Remus  600  Survey  line  of  AUVs  that  the  Navy  uses  extensively  have  an 
ll.SkWh  battery  paek  whieh  allows  for  mission  times  of  at  least  33  hours  [8].  At  that 
power  rate,  the  MBARI  buoy  eould  provide  the  power  needs  of  one  of  these  AUVs 
indefinitely  in  normal  operating  conditions,  and  if  a  similar  setup  were  employed  by  the 
Navy  it  has  the  potential  to  eliminate  the  AUVs  dependence  on  a  support  ship.  The 
drawback  is  that  a  bright  yellow  buoy  is  not  practical  in  every  situation  where  the  Navy 
may  want  to  deploy  one  of  their  AUVs.  Other  designs  of  wave  energy  converters  lend 
themselves  better  to  being  deployed  in  a  non-obvious  fashion  if  that  is  important;  see 
Chapter  II  for  more  details  on  wave  energy  converter  designs. 

A  persistent  presence  of  these  types  of  devices  without  a  need  for  constant 
support  vessel  use  is  a  game  changing  idea.  Whether  providing  surveillance  for  the 
defense  of  U.S.  ports  and  bodies  of  water  or  gathering  intelligence  in  areas  of  interest,  the 
benefits  gained  through  the  use  of  minimally  supported  AUVs  would  be  significant. 

C.  THESIS  OBJECTIVES,  CONTRIBUTIONS  AND  ORGANIZATION 

The  expressed  interest  of  this  research  is  to  aid  in  the  understanding  of  linear 
generators  for  use  by  the  DOD.  In  particular,  what  are  the  benefits  and  limitations  of 
linear  generators  and  where  is  it  appropriate  to  commit  further  research  for  their  use.  In 
order  to  accomplish  this  goal  a  Simulink  model  is  developed  that  allows  for  a  flexible 
study  of  linear  generator  behavior  in  different  configurations.  The  model  is  used  to 
explore  the  feasibility  of  differing  input  functions  as  well  as  the  effects  of  changing 
certain  parameters  of  a  model  generator. 

This  work  is  separated  into  five  chapters.  Chapter  I  discusses  background 
information  that  is  germane  to  the  research  in  order  to  aid  the  reader  in  understanding  the 
motivations  and  implications  of  the  work.  Chapter  II  is  a  dialog  on  technologies  that  are 
related  to  the  real  world  application  of  linear  generators.  Chapter  III  speaks  to  the  methods 
and  specific  scenarios  used  in  developing  various  linear  generator  models.  In  Chapter  IV, 
the  results  from  the  simulations  are  discussed  and  the  meanings  behind  said  results  are 
teased  out.  Chapter  VI  draws  conclusions  on  the  work  and  discusses  the  need  for  future 
work.  All  MATLAB  and  Simulink  code  is  included  in  its  entirety  in  the  Appendix. 
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II.  RELATED  TECHNOLOGY:  POWER  ELECTRONICS  AND 

ENERGY  HARVESTERS 


This  chapter  will  discuss  technologies  that  are  related  to  power  generation 
utilizing  linear  generators.  This  ineludes  a  look  at  power  eleetronies  sueh  as  solid  state 
power  eonverters  and  battery  management  systems.  There  will  also  be  a  diseussion  of 
emerging  methods  of  power  generation  that  are  applieable  to  linear  generators  sueh  as 
energy  harvesting  of  meehanieal  vibrations  and  various  wave  energy  eonverters. 

A,  SOLID  STATE  POWER  CONVERTERS 

The  reeent  progress  in  brushless  de  maehines  would  not  be  possible  without  the 
use  of  advaneed  solid  state  power  eonverters.  In  order  to  allow  for  a  dc  maehine  to 
operate  without  brushes,  the  mechanical  commutation  that  was  being  performed  by  the 
brushes  must  be  replaeed  with  eleetronies  that  perform  the  same  function.  When  this  is 
done  the  maehine  is  said  to  be  eleetrically  eommutated. 

Solid  state  power  eonverters  exist  in  a  large  variety  of  topologies.  Generieally, 
they  fall  into  four  groups:  ac-ae  eonverters,  de-de  eonverters,  ae-dc  converters  and  de-ac 
eonverters.  It  is  also  possible  for  a  single  eonverter  to  be  designed  to  operate  in  both 
direetions.  Within  eaeh  of  these  groups  the  power  supplies  ean  be  again  split  into  two 
additional  types:  linear  and  switehing  power  supplies. 

Linear  power  supplies  operate  the  transistors  in  linear  mode  in  order  to  provide 
the  desired  output  power  level  and  quality.  The  advantage  of  linear  power  supplies  is 
their  performanee.  The  disadvantage  of  linear  power  supplies  is  their  generation  of  exeess 
heat  and  inherent  need  to  dissipate  it,  this  heat  is  also  a  form  of  wasted  power  and  thus, 
linear  power  supplies  have  lower  peak  effieieneies  than  switehing  power  supplies  [10]. 

Switehing  power  supplies  operate  by  turning  a  solid  state  switeh,  whieh  eonneets 
the  souree  power  to  the  output  power,  on  and  off  at  a  high  frequency.  By  using  an 
inductance  (L)  or  induotanee-oapaoitanee  {LC)  cireuit  design  the  output  voltage  ean  made 
lower  or  higher  than  the  souree  voltage.  The  disadvantages  of  switehing  power  supplies 
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are  added  eomplexity,  whieh  ean  lead  to  higher  eost  and  lower  reliability.  The  advantages 
to  switching  power  supplies  are  their  efficiency,  weight  and  size. 

By  regulating  their  duty  cycle,  switching  power  supplies  can  continuously  correct 
to  provide  a  regulated  output  voltage  from  a  non-regulated  source.  In  this  way,  one  could 
charge  a  battery  at  a  precise  voltage  or  a  precise  current,  even  using  a  power  source 
which  fluctuates.  The  reverse  is  also  true;  switching  power  supplies  can  provide  a 
variable  output  utilizing  a  constant  power  source. 

A  power  supply  that  takes  a  variable  input  voltage  (usually,  but  not  necessarily 
periodic)  and  creates  a  dc  output  is  called  a  rectifier.  Most  of  the  power  supplies  that  are 
associated  with  small  portable  electronics  are  this  type.  The  output  of  a  rectifier  can  be 
regulated  or  unregulated,  depending  on  the  quality  of  the  source  and  the  quality  desired  in 
the  output.  For  sensitive  equipment,  such  as  computers,  the  output  will  be  regulated.  In 
less  sensitive  applications,  like  simple  battery  chargers,  the  output  may  be  less  regulated 
to  reduce  complexity  and  cost. 

The  Institute  of  Electrical  and  Electronics  Engineers  (IEEE)  defines  a  power 
supply  that  takes  a  dc  input  and  creates  an  ac  output  as  an  inverter  [11].  Inverters  range 
from  simple  oscillators  to  sophisticated  solid  state  electronic  circuits.  An  inverter  does 
not  have  to  create  a  sine  wave  output,  it  could  create  a  triangle  wave,  square  wave  or  any 
other  periodic  output  as  this  would  still  meet  the  definition  of  an  ac  output.  As  with  the 
other  power  supplies,  the  degree  of  complexity  depends  on  the  application.  A  common 
means  of  quantifying  the  quality  of  an  inverter  that  produces  a  sine  wave  output  is  to 
measure  its  total  harmonic  distortion  (THD),  where  THD  is  a  ratio  of  the  root  mean 
square  (rms)  of  the  amplitude  of  all  of  the  distortion  harmonics  over  the  rms  of  the 
amplitude  of  the  fundamental  [10].  Companies  often  specify  requirements  for,  or  the 
quality  of,  their  products  in  this  manner. 

B,  BATTERY  CHARGERS  AND  MANAGEMENT  SYSTEMS 

The  variety  of  batteries  available  on  the  market  is  seemingly  endless,  and  while 
not  all  batteries  need  a  charger  designed  specific  to  their  specifications,  quite  a  large 
number  of  them  do.  There  are  then,  appropriately,  a  large  number  of  battery  chargers 
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available  as  well.  Despite  these  large  numbers,  these  chargers  can  be  reduced  to  two 
categories:  simple  (dumb)  chargers  and  smart  chargers. 

Simple  chargers  supply  a  constant  output  voltage  without  regard  to  the  state  of  the 
battery.  This  output  can  be  a  constant  dc  or  a  pulsed  output  that  remains  at  a  set  value. 
Because  simple  chargers  lack  any  ability  to  sense  the  condition  of  the  battery,  they 
typically  charge  at  a  slower  rate  (utilizing  a  lower  voltage)  to  prevent  severe 
overcharging  which  can  permanently  damage  a  battery.  Simple  chargers  are  generally 
cheaper  which  makes  them  an  attractive  option  under  certain  circumstances  [12]. 

Smart  chargers  vary  their  output  based  on  the  conditions  of  the  battery.  Some 
smart  chargers  are  designed  to  be  used  with  smart  batteries  which  have  onboard  sensors 
which  communicate  to  the  charger,  while  other  smart  chargers  operate  without  the  luxury 
of  on-battery  sensors.  For  those  operating  without  the  use  of  on-battery  sensors,  they 
monitor  their  own  output  voltage  and  current  to  determine  the  state  of  the  battery,  and  in 
many  cases  this  method  is  sufficient  and  safe.  In  the  case  of  lithium  batteries;  however,  it 
is  insufficient,  which  is  why  lithium  ion  batteries  always  contain  onboard  sensors  in  order 
to  communicate  their  state  to  their  charger  [12]. 

When  a  smart  charger  acts  to  maintain  a  battery,  whether  it  is  a  single  cell  or  a 
pack,  within  a  set  of  parameters  it  is  performing  the  duties  of  a  battery  management 
system  (BMS).  This  is  necessary  because  many  modem  battery  types  become  unstable 
and  seriously  dangerous  outside  of  a  small  window  of  conditions.  According  to  [12]  a 
BMS  may  perform  any  of  the  following  functions: 

•  Monitor  the  battery; 

•  Protect  the  battery; 

•  Estimate  the  battery’s  state; 

•  Maximize  the  battery’s  performance; 

•  Report  status  and/or  history  of  the  battery  to  users  or  devices. 

Some  devices  have  only  single  cell  batteries  (like  cell  phones),  while  other 
devices  utilize  a  battery  pack  which  contains  multiple  cells  arranged  in  series  or  parallel 
or  a  combination  of  the  two.  Where  the  technology  of  BMSs  is  concerned,  there  is  a 

sizeable  leap  in  complexity  associated  with  charging  a  battery  pack  vice  a  single  cell. 
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This  is  because  individual  cells  still  need  individual  attention  in  order  to  charge  quickly 
and  safely  and  maintain  their  health.  As  a  result,  these  BMSs  may  have  entire  arrays  of 
sensors  and  may  require  multiple  outputs  to  be  tailored  to  different  groups  or  individual 
cells. 

Regulated  power  supplies  enable  complex  management  systems  created  for  the 
purpose  of  safely  and  efficiently  charging  batteries.  It  is  this  pair  of  technologies  that 
could  enable  an  inexpensive  linear  generator  with  low  power  quality  output  to  be  used  to 
supply  power  to  a  bus  or  a  battery  bank. 

C.  APPLICATIONS  OF  LINEAR  GENERATORS 

Linear  generators,  also  called  linear  alternators,  have  a  growing  list  of  modem 
applications.  Linear  generators  are  particularly  well  suited  to  cases  which  involve 
reciprocating  motion.  Some  of  the  most  prominent  or  promising  applications  are  in  firee- 
piston  Stirling  engines,  crankless  internal  combustion  engines  and  energy  harvesting, 
including  wave  energy  converters  and  vibration  harvesting  attenuators. 

1.  Free-piston  Stirling  Engine 

Free-piston  Stirling  engines  have  become  an  emerging  technology  use  of  the 
linear  generator.  Stirling  engines  require  only  a  difference  in  temperature  to  operate, 
which  allows  them  to  be  useful  in  a  vast  number  of  situations.  They  can  operate  with  a 
traditional  fuel  source  (i.e.,  burning  of  fossil  fuels)  or  more  extraordinary  heat  sources 
like  radioisotope  decay  or  concentrated  sunlight;  they  are  also  useful  for  utilizing  waste 
heat  from  various  processes  [13]. 

Two  uses  stand  out  in  the  realm  of  free-piston  Stirling  engines.  The  first  is  a  space 
application,  the  Stirling  radioisotope  generator  (SRG),  shown  in  Figure  2.  The  SRG  is  a 
new  and  more  efficient  take  on  the  radioisotope  thermoelectric  generator  (RTG)  that 
NASA  has  been  using  for  years.  While  the  RTG  achieves  an  efficiency  of  around  six 
percent,  the  new  SRG  will  be  four  times  as  efficient  and  therefore  use  only  a  fourth  of  the 
plutonium  required  of  the  older  units  [14]. 
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Figure  2.  Free-piston  Stirling  engine  and  linear  generator  portions  of  an  SRG, 

from  [14], 


The  seeond  applieation  of  the  free-piston  Stirling  engine  is  its  employment  in 
solar  eoneentrators.  When  utilized  with  solar  concentrators,  the  power  of  the  sun  is 
reflected  by  a  parabolic  dish  or  trough  and  focused  on  a  free-piston  Stirling  engine,  which 
operates  a  linear  generator  to  produce  electricity  from  solar  energy  without  the  need  of 
expensive  photovoltaic  cells;  a  depiction  of  this  setup  can  be  seen  in  Figure  3.  In  this 
configuration,  solar  energy  can  be  collected  at  an  efficiency  of  nearly  30  percent  [15]. 


Figure  3.  Depiction  of  a  solar  concentrator  utilizing  a  reflective  parabolic  dish 

and  a  Stirling  engine,  from  [15]. 
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2,  Crankless  Internal  Combustion  Engines 

In  this  clever  applieation  of  the  linear  generator  the  standard  internal  eombustion 
engine  is  reimagined  to  replaee  the  pistons  with  linear  generators.  The  resulting  maehine 
removes  the  need  for  gear  boxes,  drive  shafts  and  additional  generators.  The  benefits  are 
three  fold.  The  first  benefit  is  a  reduetion  of  size  and  weight  whieh  adds  to  the  overall 
effieieney  and  responsiveness  of  a  vehiele.  The  seeond  benefit  is  a  reduetion  in  moving 
parts  whieh  ean  lead  to  lower  eost,  greater  simplieity  and  higher  reliability  in  basie 
applieations.  The  third  benefit  is  an  inerease  in  effieieney  by  eliminating  unneeessary 
mechanieal  losses  eaused  by  the  outmoded  mechanieal  systems.  An  example  free-piston 
linear  generator  designed  for  a  hybrid  vehiele  applieation  ean  be  seen  in  Figure  4,  where 
TDC  stands  for  top  dead  eenter  and  BDC  for  bottom  dead  eenter,  representing  the 
respeetive  extremes  of  the  piston  stroke.  The  spring  baek  volume  aets  as  an  air  spring  that 
limits  the  motion  of  the  piston  on  one  end  of  the  eycle. 


Figure  4.  Example  of  a  hybrid  vehicle  free-piston  linear  generator,  from  [16]. 


This  application  of  the  linear  generator  is  espeeially  relevant  given  the  growing 
trend  of  vehicle  manufactures  to  incorporate  eleetric  drive  technologies  into  their 
vehieles.  Vehieles  sueh  as  the  BMW  i3  currently  utilize  range  extenders  with  no  direet 
eonneetion  to  the  drivetrain,  instead  these  internal  eombustion  engines  solely  power 
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generators  whieh  in  turn  provide  power  solely  to  eliarge  the  vehiele  battery  and  the 
eleetrie  drive  motors  [17]. 

3,  Vibration  Harvesters 

Vibration  harvesters  are  a  niehe  market  of  speeialized  deviees  whieh  generate 
eleetrie  energy  from  wasted  meehanieal  energy  in  the  form  of  vibration.  Most  eommonly, 
these  deviees  utilize  piezoeleetrie  material  to  eonvert  the  mechanical  vibration  into 
electric  energy  utilizing  the  piezoelectric  effect,  whereby  an  applied  mechanical  stress 
generates  an  electric  charge.  This  very  low  power  method  of  energy  harvesting  has 
applications  mostly  in  remote  sensors  where  electrical  power  is  not  available  (e.g., 
wireless  sensors  on  train  cars).  While  piezoelectric  harvesters  are  primarily  utilized  in 
situations  with  high  frequency  vibrations,  situations  with  relatively  low  frequency 
vibrations  lend  themselves  better  to  small  linear  generators.  Companies  like  LORD 
MicroStrain  produce  both  types  of  technologies  commercially  for  use  in  their 
accompanying  sensor  systems  [18],  [19],  [20]. 

While  vibration  harvesting  can  be  very  useful  for  providing  power  to  remote 
sensors,  using  a  linear  generator  to  perform  this  action,  when  appropriate,  can  realize 
several  advantages.  Linear  alternators  have  the  potential  to  provide  more  power  than 
piezoelectric  systems  given  a  large  enough  mechanical  input.  More  interestingly,  though, 
linear  generators  can  provide  controlled  mechanical  damping  in  a  system.  Many  existing 
damper  applications  could  benefit  from  replacing  existing  dampers  with  computer 
controlled  active  dampers,  such  as  replacing  the  shock  absorbers  in  a  vehicle  with 
compact  linear  generators.  A  U.S.  patent  was  awarded  for  just  such  an  application  in 
October  of  2005  [21]. 

4,  Wave  Energy  Converters 

A  wave  energy  converter  (WEC)  is  a  device  which  generates  electrical  power  by 
capturing  the  power  found  in  ocean  waves.  WECs  are  grouped  into  several  distinct 
categories: 

1 .  Surface  point  absorbers — are  typically  buoys  which  harness  power  from 

the  up-down  motion  they  see  while  riding  atop  waves. 
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2.  Surface  attenuators — are  floating  snake-like  deviees  whieh  harness  power 
through  flexing  aetion  as  waves  travel  along  their  length. 

3.  Subsurfaee  point  absorbers — aet  like  surfaee  point  absorbers  exeept  that 
they  follow  the  rise  and  fall  of  the  oeean  surfaee  not  through  buoyaney, 
but  through  the  differenee  in  pressure  that  is  eaused  by  the  passing  waves. 

4.  Overtopping  deviees — are  large  deviees  whieh  allow  incoming  waves  to 
overtop  them  and  eause  the  water  height  in  them  to  exeeed  the  water 
height  of  the  water’s  surfaee.  This  resulting  potential  energy  is  harnessed 
by  allowing  the  higher  level  water  to  drain  through  turbines. 

5.  Wave  surge  eonverters — are  deviees  whieh  harness  the  side-to-side 
oseillations  that  oeeur  in  shallow  water  waves.  These  deviees  are  often 
large  flaps  that  are  hinged  at  the  sea  floor. 

Of  all  of  these  eategories  of  deviees,  only  overtopping  deviees  do  not  lend 
themselves  to  the  direet  applieation  of  linear  generators.  It  may  be  surprising  then  to 
know  that  the  majority  of  existing  wave  energy  eonverter  projeets  utilize  hydraulies  to 
store  power  whieh  they  then  use  to  spin  turbines  attached  to  generators.  This  is  partially 
beeause  the  higher  the  foree  and  the  lower  the  speed,  the  more  suited  the  applieation  is  to 
a  hydraulie  power  take-off  method  [22].  As  the  teehnology  matures,  however,  there  is  a 
developing  affinity  for  direet  drive  linear  generators  as  the  power  take-off  method  of 
ehoiee  in  wave  energy  eonverters.  This  growing  fondness  for  linear  generators  is 
espeeially  true  in  smaller  deviees  in  the  kilowatt  range,  and  they  are  even  being 
implemented  in  mueh  larger  deviees,  sueh  as  the  2MW  Arehimedes  Wave  Swing  shown 
in  Figure  5. 

There  are  ehallenges  still  to  be  solved  in  the  design  and  applieation  of  linear 
generators  on  the  seale  required  for  eommereial  energy  production.  But  the  allure  of 
lower  maintenanee  requirements  and  higher  reliability  on  objeets  that  reside  below  the 
oeean  surfaee  is  sure  to  keep  linear  generators  a  topie  of  interest  and  development  in  the 
wave  energy  converter  eommunity. 
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Figure  5.  Permanent  magnet  linear  generator,  with  workers  inside,  found  in  the 

Archimedes  Wave  Swing,  from  [22], 
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III.  LINEAR  GENERATOR  SYSTEM  DESIGN  AND  MODELING 


Chapter  III  will  explain  some  basic  machine  theory  necessary  for  the 
understanding  of  the  constructed  model.  Next  linear  synchronous  machines  will  be 
covered  and  finally  the  chapter  will  go  into  detail  regarding  construction  of  the  Simulink 
model.  Some  assumptions  are  made  about  the  construction  of  the  machine  in  order  to 
make  the  model  more  manageable.  These  lead  to  some  limitations  in  its  use,  however, 
and  present  room  for  further  work. 

Although  non-sinusoidal,  the  model  machine  is  assumed  to  be  half-wave 
symmetric  which  prevents  it  from  containing  any  even  harmonics.  It  is  further  assumed 
that  the  model  generator  does  not  enter  saturation  at  any  point  in  the  simulation  and  that  it 
does  not  exhibit  any  demagnetization  by  the  permanent  magnets.  Cogging  torque  has  also 
been  neglected,  which  may  be  of  some  significance  in  a  linear  machine  that  is  operating 
at  low  frequencies. 

Sign  and  naming  conventions  for  the  model,  as  well  as  the  detailed  governing 
equations,  are  borrowed  largely  from  Chapmen  et  al  [23].  Attempts  are  made  to  point  out 
where  these  conventions  may  differ  from  the  standard  or  where  there  may  not  be  a 
standard  at  all. 

A,  SYNCHRONOUS  MACHINES 

A  synchronous  machine  is  one  in  which  the  physical  speed  of  the  machine  is  an 
exact  multiple  of  the  electrical  frequency  of  the  machine.  In  these  machines  the  magnetic 
fields  of  the  stator  (the  stationary  part  of  the  machine)  and  rotor  (the  rotating  part  of  the 
machine)  rotate  at  the  same  speed  and  stay  aligned  with  each  other  under  normal 
operation.  While  it  is  not  fixed,  and  either  the  rotor  or  stator  may  occupy  the  inner  portion 
of  the  machine,  it  is  more  common  to  see  the  rotor  on  the  inside  surrounded  by  the  stator. 
The  magnetic  poles  produced  by  the  electrical  windings  are  perpendicular  to  the 
windings.  Often  in  figures,  such  as  Figure  6,  the  machine  windings  are  depicted  as  a 
single  lump  winding.  In  modem  machines,  however,  this  is  almost  never  the  case; 
instead,  distributed  windings  are  employed,  with  the  distribution  done  in  a  manner  to 
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produce  the  baek  eleetro motive  foree  (EMF)  shape  desired  (e.g.,  sinusoidal  back  EMF  or 
trapezoidal  baek  EMF). 

When  there  is  no  angle  between  the  magnetie  poles  in  the  rotor  and  stator  there  is 
no  torque  produced,  but  as  the  angle  gets  larger  it  produces  greater  torque  until  it  reaehes 
90°  where  it  produees  the  maximum  torque.  Adding  any  load  beyond  this  maximum 
torque  will  eause  the  maehine  to  produee  less  torque  and  will  induee  an  unstable  behavior 
where  the  torque  angle  will  quiekly  grow.  Beyond  180°,  the  torque  flips;  if  this  oecurs  in 
an  eleetrieal  maehine,  the  maehine  is  said  to  have  stalled.  Conditions  of  maximum  and 
minimum  torque  are  depleted  in  Figure  7.  Also  shown  in  Figure  6  is  the  rotor  angle  Or.  It 
should  be  noted  that  in  a  machine  with  more  than  two  poles  instability  will  be  reaehed  at 
angles  less  than  90°. 


Figure  6.  Simple  two  pole  eleetrieal  maehine  diagram  showing  the  rotor  and 
stator,  (a)  End  view;  (b)  eutaway  eross-seetion  view,  from  [24]. 


Figure  7.  Magnets  produeing  (a)  maximum  torque  and  (b)  zero  torque. 
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In  terms  of  a  synchronous  motor,  the  rotating  magnetic  field  of  the  stator  is 
generated  by  a  periodic  current  waveform  supplied  from  a  power  source.  The  magnetic 
field  strength  //  in  a  wire  in  free  space  is 


Ijir 


(1) 


where  /  is  the  current  in  the  wire  and  r  is  the  radius  from  the  center  of  the  wire.  The  flux 
density  B  is  related  to  H  by 

B  =  jliH.  (2) 


Magnetic  permeability  //  is  the  measure  of  a  material’s  ability  to  support  a  magnetic  field. 
Normally,  //  is  expressed  as  the  product  of  the  permeability  of  free  space  constant  jig  (e.g., 
the  permeability  of  a  vacuum)  and  the  dimensionless  relative  permeability  of  the 
substance  jir,  ■  [24] 

The  magnetic  field  in  the  rotor  of  a  synchronous  machine  can  either  be  produced 
the  same  way,  using  an  electrical  current,  or  by  placing  permanent  magnets  on  the  rotor, 
thus  eliminating  the  need  for  electrical  connections  to  the  rotating  portion  of  the  machine. 
In  terms  of  a  generator,  it  is  the  rotor  that  produces  the  rotating  magnetic  field,  either 
through  use  of  electrical  current  or  permanent  magnets.  Each  winding  in  the  stator  then 
sees  a  continuously  varying  magnetic  field,  which  by  the  laws  of  magnetic  induction 
induce  a  voltage  on  the  winding.  In  1833,  Heinrich  Lenz  first  stated  that  motors  and 
generators  were  reversible;  his  findings  were  later  summarized  by  James  Clerk  Maxwell 
in  Lenz’s  law,  which  states  that 


dl 

s  = - 

dt 


(3) 


where  e  is  the  EMF  measured  in  volts,  X  is  the  flux  linkage,  d  is  the  derivative  operator 
and  t  is  time.  The  flux  linkage  of  a  winding  is  related  to  the  vector  of  the  magnetic  field 
strength  B  by 

/I  =  j  B.  (4) 
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where  •  is  the  dot  produet  operator  and  variables  in  bold  represent  veetors.  The 
integration  is  taken  over  the  area  inside  the  winding.  This  ean  be  simplified  in  the  ease  of 
eoils  of  wire  to 

X^BA  (5) 

where  A  is  the  surfaee  area  inside  the  eoil.  If  there  are  N  turns  of  wire  in  the  same  area 
with  the  same  B  field  aeting  on  them,  then  the  total  flux  linkage  is 

A  =  NBA.  (6) 

To  find  the  total  voltage  produced  in  a  stator  winding  we  need  to  look  at  both  the 
EMF  and  the  electrostatic  voltage  in  the  winding.  The  total  voltage  in  the  winding  then  is 


dX 


(V) 


where  r  and  i  are  the  resistance  and  current  in  the  stator.  In  a  three-phase  machine  with 
stator  phases  a,  b  and  c,  the  voltage  in  the  stator  can  be  described  by 


'  abcs 


r\  +  — k 

spates  j,  ^abes^ 

at 


(8) 


where  Yabes,  iabes  and  Xabes  are  generalized  in  the  form 


^abes 


fa. 

fb. 
fes 


(9) 


and/as,  fbs  and  fc.  are  the  a,  b  and  c  phases  of  stator  voltage,  current  or  flux  linkage. 

In  a  machine  with  a  sinusoidal  back  EMF,  the  total  flux  linking  the  rotor  and 
stator  can  be  found  by  combining  the  flux  from  the  permanent  magnets  on  the  rotor  and 
the  flux  from  all  the  windings.  The  flux  experienced  at  a  given  winding  from  the 
permanent  magnet  is  dependent  on  the  rotor  position  at  that  time.  More  specifically,  it  is 
dependent  on  the  rotor  electrical  angle  Ore-  The  electrical  position  of  the  rotor  Ore  is 
related  to  the  mechanical  rotor  angle  Orm  by  the  number  of  magnetic  poles  p  in  the  motor 
per  phase,  where 
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(10) 


e,. 


e,. 


p. 

2  ■ 


The  flux  linkage  matrix  "kabcs  then  ean  be  defined  by 


abcs 


(11) 


where  the  peak  flux  linkage  due  to  the  permanent  magnet  is  a  eonstant  sinee  it  is 

assumed  that  the  permanent  magnets  are  not  undergoing  any  demagnetization.  The  stator 
self-inductanee  matrix  for  a  non-salient  machine  is  given  by 


(12) 


where  Lis  and  L^s  are  the  stator  leakage  and  mutual  inductances  respectively. 

Defined  for  generator  action,  the  mechanical  input  torque  im  and  electrical 
counter-torque  le  can  be  related  to  mechanical  rotor  speed  comi  by 


T  —T  =J  —  O)  +B  CO  , 

e  m  i_,  rm  m  rm‘ 

at 


(13) 


where  J  is  the  rotational  inertia  of  the  rotor  and  is  the  mechanical  damping  due  to 
friction.  The  Jterm  takes  into  account  the  machine’s  resistance  to  a  change  in  speed  and 
the  Bm  term  accounts  for  losses  in  the  machine  due  to  friction,  both  in  bearings  and  due  to 
air  resistance. 


The  electric  torque  Xe  that  a  machine  generates  is  specific  to  the  geometry  of  that 
machine.  It  can  be  complicated  by  a  machine  that  has  a  non-sinusoidal  winding 
distribution  or  by  a  machine  that  is  salient  (the  width  of  the  air  gap  between  the  stator  and 
rotor  varies  because  one  or  both  is  not  round).  If  we  assume  the  machine  is  sinusoidal  in 


21 


its  winding  distribution  and  has  a  constant,  non-salient,  air  gap  then  the  eleetrieal  torque 
equation  can  be  expressed  by 


The  eleetrieal  machine  that  is  depicted  in  the  model  will  be  of  a  non-sinusoidal 
nature  and  the  equations  specific  to  that  problem  will  be  discussed  further  in  the 
modeling  section. 

B.  LINEAR  SYNCHRONOUS  MACHINES 

A  linear  machine  is  coneeptually  the  same  as  a  traditional  rotating  machine.  It  ean 
be  described  as  eutting  and  unrolling  a  rotating  machine  so  that  it  is  flat.  A  rotating 
machine  is  depicted  in  Figure  8,  with  the  same  machine  unrolled  in  Figure  9  where 

stator  windings  depicting  current  flow  into  the  page; 
a[,  b[,  c'  stator  windings  depicting  current  flow  out  of  the  page; 

rotor  windings  depicting  current  flow  into  the  page; 
rotor  windings  depicting  current  flow  out  of  the  page; 

(p^  angular  displacement  along  the  rotor  eireumferenee; 

(p^  angular  displaeement  along  the  stator  circumference; 

co^  rotor  speed. 
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Figure  8.  Example  two-pole,  three-phase  salient-pole  synchronous  machine 
showing  the  general  layout  of  windings  and  major  axis,  from  [23]. 


axis  cs  axis  bs  axis  as  axis 


Figure  9.  Rotating  to  linear  machine  unrolling  concept  visualized.  The 
fluctuating  line  represents  the  salient  air  gap,  from  [23]. 
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It  is  common  to  model  or  discuss  a  linear  maehine  using  the  same  equations  and 
variables  as  a  rotating  machine.  The  results  must  be  translated  to  terms  and  variables  that 
ean  represent  a  linear  machine.  Torque  can  be  related  to  foree  by 

P^tco^^^FV,  (15) 

where  P  is  power,  F  is  force  and  V  is  linear  veloeity.  Electrieal  power  into  the  three  phase 
machine  P  can  be  ealculated  using  machine  parameters  by 

P^^asKs+^bshs+^cshs-  (16) 


The  length  of  the  linear  machine  //„  can  be  approximated  using  the  radius  of  the 
rotating  machine  by  //mS27iR.  This  is  not  exact  because  one  has  to  choose  a  radius,  as  the 
radius  of  the  rotor  and  stator  will  differ.  In  the  case  of  small  maehines  the  difference 
could  amount  to  a  large  pereentage  of  the  total  circumference.  Unlike  a  rotating  machine 
where  the  entirety  of  the  stator  always  is  interaeting  with  the  rotor  and  likewise  the 
entirety  of  the  rotor  is  always  interacting  with  the  stator,  portions  of  a  linear  machine  are 
not  interaeting  with  each  other  at  any  given  time.  This  is  because  either  the  stationary 
portion  of  the  maehine  (still  referred  to  as  a  stator)  or  the  moving  portion  of  the  machine 
(referred  to  as  a  translator)  is  designed  longer  than  the  other.  This  extra  length  is  the 
allowable  travel  of  the  linear  maehine  without  produeing  generally  undesirable  results, 
and  is  shown  in  Figure  10.  There  is  no  standard  for  whether  the  stator  or  the  translator  is 
the  portion  with  permanent  magnets,  the  ehoice  is  a  balanee  of  cost  and  performance. 
Permanent  magnet  material  is  more  expensive  than  wires,  but  magnets  do  not  require  any 
electrieal  connections  allowing  for  cheaper  more  reliable  designs. 


Translator  motion 


I  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  II  N  II  S  I 

o'  b  c'  a  b  c 

I - - 1 —  Overhang  -  usable  travel — | 

Figure  10.  Example  three-phase  linear  machine  setup  with  permanent  magnet 

translator. 
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c. 


MODELING 


In  [25],  Chapman  et  al.  discuss  modeling  a  three-phase  non-sinusoidal  permanent 
magnet  synehronous  maehine  using  a  rotor  reference  frame  method.  Using  the  Chapman 
and  company  machine  description  and  experimental  results,  the  following  Simulink 
model  is  developed  and  applied  to  a  linear  maehine  of  notionally  similar  eharacteristics. 

In  the  synehronous  machines  section,  machine  equations  were  presented  for, 
arguably,  the  simplest  type  of  rotating  synehronous  machine,  a  three  phase  permanent 
magnet  maehine  with  a  constant,  non-salient  air  gap  which  has  a  sinusoidal  winding 
distribution  resulting  in  a  sinusoidal  back  EMF.  The  maehine  deseribed  and  modeled  here 
(and  in  [25])  however  has  a  non-sinusoidal  back  EMF.  A  trapezoid  is  a  common 
waveform  for  a  non-sinusoidal  machine  to  utilize.  There  are  several  reasons  for  choosing 
this  including  easier,  eheaper  construetion,  inereased  power  density  and  smaller  inverter 
sizes  for  motors  [25].  The  non-sinusoidal  ehoice  results  in  machine  equations  which  must 
account  for  the  harmonic  components  of  the  waveforms.  While  the  equations  account  for 
an  infinite  number  of  harmonic  components,  in  practice  only  a  few  are  used  in  modeling. 
In  this  ease  it  is  also  assumed  that,  although  non-sinusoidal,  the  baek  EMF  is  half-wave 
symmetrie  whieh  eliminates  the  even  harmonies. 

1,  Machine  Variable  Equations  for  a  Non-Sinusoidal  Machine 

In  the  non-sinusoidal  maehine  equations  the  stator  voltage  equation  (8)  and  the 
self-inductanee  matrix  (12)  remain  the  same  while  the  flux  linkage  equation  (11) 
becomes 


'^abcs  ^i^abcs'^ 
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where  the  Kn  coefficients  from  (17)  represent  the  magnitude  of  the  harmonic  components 
of  the  flux  linkage  and  the  coefficients  K[  -  nK^  and  represent  the  harmonic  components 
of  the  torque  in  (18)  [25]. 


2,  Rotor  Reference  Frame  Model  Using  qdO  Variable 

The  Simulink  model  also  utilizes  the  rotor  referenee  frame  through  a 
transformation  of  variables.  Using  a  transformation  of  variables  in  order  to  analyze  three 
phase  systems  dates  baek  to  1929  when  Robert  H.  Park  first  published  the  Park’s 
transformation  [23].  Sinee  then,  the  very  similar  direet-quadrature-zero  (qdO) 
transformation  has  replaeed  the  Park’s  transformation  as  a  standard  transformation  for 
three  phase  systems.  In  the  qdO  transformation,  the  a,  b  and  c  phase  variables  are 
projected  onto  a  new  set  of  referenee  frame  axes  that  are  defined  as  the  q,  d  and  0  axes, 
as  shown  in  Figure  1 1 .  Where  fas,fbs  and  fcs  are  the  a,  b  and  c  axes  and  fds  and  fqs  are  the  d 
and  q  axis.  The  0-axis  is  not  shown,  but  would  originate  from  the  interseetion  of  the  other 
axes  and  be  positive  out  of  the  page. 

The  reason  the  0-axis  is  not  shown  is  because  when  the  qdO  transformation  is 
applied  to  a  three-phase  system,  whieh  has  sinusoidal  geometry  (whieh  is  eommon),  the 
0-axis  variables  go  to  zero  and  only  the  q  and  d  axes  variables  are  left.  A  major  advantage 
to  using  the  qdO  transformation  is  that  it  allows  for  simplified  direet  torque  eontrol  in 
eleetrieal  maehines. 

Despite  no  need  for  a  transformation  (as  there  is  no  controller  and  the  model  is 
non-sinusoidal),  the  model  needs  to  keep  this  format  for  added  utility  in  any  follow-on 
work  where  a  eontroller  might  be  added. 
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Figure  11.  Axis  relationships  for  an  abe  to  qdO  variable  transformation.  The  0- 
axis,  not  depicted,  originates  from  the  intersection  of  the  axes  and  is 
positive  out  of  the  page,  from  [23]. 

In  order  to  perform  a  qdO  transformation,  multiply  the  machine  variables  by  the 
transformation  matrix  K.  Since  we  are  interested  in  transforming  stator  variables  onto  the 
rotor  reference  frame,  we  will  use  an  electrical  rotor  angle  {Ore)  as  the  reference  in  the 
transformation  and  the  resulting  transformation  matrix  will  be  termed  where 


the  new  stator  variables  in  the  qdO  rotor  reference  frame  are  now  generalized  as  , 
where 

/;/ 

l:  (20) 


and 
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=  K  f 

5  C 

The  new  q  and  d  axes  in  the  rotor  reference  frame  are  shown  in  Figure  12.  It  can  be  seen 
that  the  ^-axis  and  the  J-axis  are  in  quadrature  (i.e.,  90°  apart).  Note  too  that  the  ^-axis  is 
the  one  that  is  referenced  by  6rm,  this  is  contrary  to  most  literature  where  the  J-axis  is  the 
reference.  This  has  the  effect  of  reversing  the  sine  and  cosine  terms  in  several  equations 
which  might  otherwise  be  easily  recognizable  to  those  familiar  with  electrical  machine 
analysis. 


Figure  12.  Cross  section  of  a  three-phase  non-salient  permanent  magnet 

synchronous  machine  with  the  q-axis  and  d-axis  shown,  from  [25]. 


To  convert  back  from  qdO  variables  to  abc  variables,  the  inverse  transformation 
matrix  (k^)  *  is  applied,  where 


(22) 


(23) 
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From  [25],  if  the  qdO  transformation  is  applied  to  the  maehine  equations  for  voltage,  flux 
linkage  and  torque  in  a  non-sinusoidal  maehine  then  (8)  for  voltage  becomes 


r  -r  ,  0  ^  0 

V  “  ^  ’ 

(24) 

r  ‘f  1  /■  d  n  f 

(25) 

r  .r  d  „  « 

dt 

(26) 

equation  (17)  for  flux  linkage  becomes 
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In  (16)  it  was  shown  that  power  can  be  calculated  using  abc  machine  variables; 
power  can  also  be  calculated  using  qdO  variables  where 


'’=^('';-;+v;4+2v;,,;,). 


(31) 


3,  Simulink  Model 

The  full  Simulink  model,  along  with  the  initial  conditions  file  and  the  post  script, 
are  included  in  the  appendix.  This  section  will  look  at  selected  portions  of  the  model  in 
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some  detail,  ineluding  an  explanation  of  where  equations  were  implemented  in  various 
subsections. 

In  order  to  provide  a  useful  reference,  the  top  level  block  diagram  of  the  model 
can  be  seen  in  Figure  13.  Flere  we  can  see  the  various  subsections  of  the  model  and  the 
signals  that  are  passed  amongst  them.  Blocks  that  are  labeled  “To  Workspace”  represent 
variables  which  are  passed  to  the  MATLAB  command  window  for  use  by  the  post-script 
file  in  creating  various  plots  presented  in  Chapter  IV.  It  is  not  uncommon  to  see  the  same 
variable  represented  differently  on  the  same  screen  or  across  different  screens  because  of 
restrictions  on  naming  and  repeating  variables  in  Simulink  windows. 


Figure  13.  Simulink  top  level  block  diagram  of  the  model. 


One  of  several  qdO  transformations  is  depicted  in  Figure  14.  This  has  been  done 
through  the  implementation  of  Equations  (19)-(21),  and  in  this  instance  the  voltage 
variables  were  transformed.  Note  that  9^^  is  used  in  this  transformation.  By  orienting  on 
the  rotor  electrical  position  in  the  transform  we  have  implemented  the  rotor  reference 
frame,  which  is  rotating  at  .  This  is  sometimes  called  the  synchronous  reference  frame 
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because  it  is  rotating  at  the  electrical  angular  velocity  of  the  air-gap  rotating  magnetic 
field.  Other  reference  frames,  such  as  the  stationary  reference  frame  or  the  mechanical 
rotor  reference  frame,  are  sometimes  used  depending  on  the  application.  Other  reference 
frames  may  rotate  at  a  speed  unrelated  to  the  machine  and  are  termed  arbitrary  reference 
frames.  Commonly,  these  the  transformations  are  depicted  by  different  superscripts  above 
the  variables.  There  is  no  established  notation  to  reference  so  readers  should  be  careful 
when  viewing  texts  on  the  subject  [23]. 
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Figure  14.  Simulink  model  subsection  depicting  a  qdO  transformation  via  the 

implementation  of  Equations  (19)-(21). 
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The  reverse  transform,  aeting  on  the  qdO  eurrent  variables  and  reforming  the  abc 
current  variables,  is  shown  in  Figure  15.  Here  we  have  implemented  Equations  (22)  and 
(23). 


Figure  15.  Simulink  model  subsection  depicting  a  reverse  qdO  transformation  via 

the  implementation  of  Equations  (22)  and  (23). 


The  bulk  of  the  permanent  magnet  synchronous  machine  is  implemented  as  a 
large  subsection  with  several  smaller  subsections  within  it,  the  larger  subsection  can  be 
seen  in  Figure  16. 
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Figure  16.  Simulink  model  subseetion  depicting  the  related  smaller  subsections 
within  the  permanent  magnet  synchronous  machine  block. 


Progressing  deeper  into  the  permanent  magnet  synchronous  machine,  the  flux 
linkage  and  current  are  determined  separately  for  each  of  the  qdO  variables.  Each  variable 
is  governed  by  different  but  related  equations.  For  the  variables,  is  discovered  by 

rearranging  the  equation  (24).  The  ^-axis  stator  current  is  found  using  the  newly 
developed  through  the  rearrangement  and  implementation  of  (27),  as  seen  in  Figure 
17  where  “lambda_m”  is  the  constant  and  the  gain  “l/(Lss)”  refers  to 
L  +(3/2)1  . 

ss  Is  ^  d  ms 
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Figure  17.  Simulink  model  subsection  depicting  the  development  of  and 

utilizing  Equations  (24)  and  (27)  respectively. 


The  J-axis  flux  linkage  and  current  variables  are  found  in  a  similar 
fashion  using  (25)  and  (28),  and  are  depicted  in  Figure  18.  While  and  utilize  (26) 
and  (29)  as  seen  in  Figure  19. 


Figure  18.  Simulink  model  subsection  depicting  the  development  of  A^  and 

utilizing  Equations  (25)  and  (28)  respectively. 
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Figure  19.  Simulink  model  subseetion  depleting  the  development  of  and 

utilizing  Equations  (26)  and  (29)  respectively. 


The  electrical  power  P  that  the  machine  generates  from  the  input  mechanical 
torque  applied  is  developed  from  qdO  variables  by  implementing  (31),  as  seen  in  Figure 
20. 


p_0s 


Figure  20.  Simulink  model  subsection  depicting  the  development  of  machine’s 

output  power  via  Equation  (31). 
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The  counter-electrical  torque  ,  that  is  produced  by  the  machine  in  opposition  to 
,  can  be  calculated  by  (30).  This  has  been  implemented  in  the  Simulink  model  of  the 
machine  in  Figure  21  where  gains  K3,  K5  and  K7  are  constants  representing  the 
magnitude  of  the  third,  fifth  and  seventh  harmonic  components  of  in  the  rotor 

reference  frame  and  gain  “-K-”  =  (3/l„',  /  2)[p  /  2). 


i_0s  Product 


Figure  21.  Simulink  model  subsection  depicting  the  calculation  of  utilizing 

Equation  (30). 

The  mechanical  forces  that  are  at  work  on  the  machine  are  represented  by  the 
force  balance  equation  (13).  Here  the  forces  of  friction,  torque  and  inertia  are  balanced  to 
govern  the  behavior  of  the  machine.  In  the  Simulink  model,  (13)  is  implemented  in  the 
force  balance  subsection  as  depicted  in  Figure  22.  In  addition  to  the  friction,  torque  and 
inertia  forces  in  (13),  a  lossless  centering  spring  has  been  simulated  out  of  necessity.  In 
all  linear  machines  the  translator  travel  is  a  design  constraint  that  must  be  accounted  for. 
This  can  be  accomplished  by  constraining  the  input  forces  or  through  the  design  of  the 
machine  itself.  In  some  instances,  hard  stops  are  placed  on  the  ends  of  the  machine  to 
prevent  over  travel.  Another  simple  solution  is  to  use  a  spring,  most  commonly  an  air 
spring,  whereby  the  translator  compresses  a  volume  of  air  as  it  travels  toward  the 
extremes  of  its  travel. 
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In  the  Simulink  model,  a  spring  is  simulated  whieh  applies  zero  force  when  the 
translator  is  centered  at  6^^  equal  to  tt,  where  6^^  has  a  total  allowable  travel  from  zero 

to  2n.  The  force  applied  by  the  spring  is  linear  and  equal  to  the  distance  from  center 
multiplied  by  a  spring  constant.  This  system  is  lossless  and  has  been  tuned  to  allow  for 
operation  at  the  target  frequency  for  the  simulation.  It  is  necessary  to  have  a  spring  since 
the  input  force  is  a  torque  that  is  independent  of  speed,  therefore  as  the  machine  speeds 
up  more  mechanical  power  is  applied  since  =  co^z^ .  If  there  is  even  the  slightest 

imbalance  where  the  machine  goes  faster  in  one  direction  that  in  the  reverse  direction,  the 
machine  becomes  unstable  and  will  exhibit  an  underdamped  response,  quickly  exceeding 
the  translator  allowable  travel  in  one  direction.  The  spring  force  is  highlighted  in  blue  in 
Figure  22. 

In  accordance  with  (13),  the  Simulink  model  is  defined  for  motor  action  and 
because  of  this,  the  summing  block  in  Figure  22  represents  as  positive  and  and  the 

mechanical  damping  as  negative.  This  frame  of  reference  is  also  taken  into 

account  when  using  the  model  to  represent  a  generator  and,  because  of  this,  is  made 

negative  as  seen  in  Figure  23.  Some  of  the  results  in  Chapter  IV  have  also  had  their  signs 
flipped  to  make  their  values  more  intuitive. 

The  output  of  the  generator  is  three  separate  phases  that  have  no  common  node. 
Because  of  the  lack  of  a  common,  the  three  phases  are  not  applied  to  a  single  three  phase 
rectifier.  Instead,  each  phase  is  applied  to  a  full  bridge  diode  rectifier  and  the  output  of 
each  rectifier  is  tied  to  a  common  dc  bus.  This  setup  is  illustrated  in  Figure  24.  The 
internals  of  the  diode  rectifiers  are  viewable  in  the  Appendix  will  and  not  be  discussed 
here. 
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Figure  22.  Simulink  model  subseetion  depicting  the  balance  of  mechanical 
forces  on  the  machine  employing  (13).  A  lossless  centering  force  is 

accented  in  blue. 


Figure  23.  Simulink  model  subsection  showing  an  example  of  the  construction 

of  for  one  setup  of  the  model. 
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Figure  24.  Simulink  model  subsection  representing  three  diode  rectifiers,  one  for 

each  phase. 
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IV.  RESULTS 


Chapter  IV  will  discuss  the  results  produced  by  the  Simulink  model  of  the 
permanent  magnet  synchronous  linear  generator,  as  well  as  look  at  some  experimental 
results  gathered  from  a  physical  brushless  dc  rotary  motor  acting  as  a  generator. 

A,  EXPERIMENTAL  RESULTS 

The  primary  purpose  of  this  section  is  to  demonstrate  a  machine,  which  is 
designed  to  function  as  a  motor,  operating  normally  as  a  generator.  For  this  experiment,  a 
delta-connected  eight-pole  permanent  magnet  brushless  dc  motor  rated  for  24V  at  4,000 
RPM  was  used.  Brushless  dc  motors  are  actually  a  misnomer;  the  motors  themselves 
cannot  operate  on  a  dc  voltage.  In  actuality,  they  are  ac  motors  which  have  separate 
power  electronics  inverters  on  the  inputs  which  generate  ac  input  signals  that  may  be 
from  a  dc  input.  Commonly,  these  motors  have  sensors  onboard  to  assist  the  power 
electronics  in  generating  the  proper  inputs  to  control  the  devices.  In  the  case  of  the 
physical  machine  used  here,  it  was  equipped  with  six  leads,  three  input-power  leads,  one 
for  each  of  the  three  phases,  and  three  sensor  leads  to  report  the  status  of  the  machine  to 
the  separate  power  electronics. 

When  operated  as  a  generator,  this  physical  machine  becomes  a  permanent 
magnet  synchronous  generator,  similar  to  the  model  generator.  The  differences  are  that 
the  physical  machine  has  sinusoidal  back  EMF  while  the  model  machine  has  non- 
sinusoidal  back  EMF  and  the  physical  machine  is  delta  connected  while  the  model 
machine  has  an  ungrounded  wye  connected  design. 

The  physical  machine  feeds  a  three-phase  full  bridge  diode  rectifier  (also  called  a 
six -pulse  rectifier)  which  is  configured  for  a  delta  connected  machine,  with  the  output  of 
the  rectifier  applied  to  a  50Q  resistor.  The  unrectified  voltage  Va  is  measured  across  the  a- 
phase  (Oa)  while  the  rectified  load  voltage  vi  is  measured  across  the  load  resistor.  Both 
measurements  are  floating  measurements  using  isolated  voltage  probes  which  do  not 
share  a  common  ground.  A  current  measurement  is  taken  after  node  A  and  measures  a 
portion  of  both  the  a  and  c-phases  The  physical  machine  setup  is  depicted  in  Figure  25. 
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Figure  25.  Three  phase  delta  connected  brushless  dc  motor  and  diode  rectifier 

experimental  setup. 


The  input  phase  voltages  and  resulting  output  voltages  for  a  theoretical  six-pulse 
rectifier  are  shown  in  Figure  26.  Part  (a)  shows  the  theoretical  three-phase  sinusoidal- 
input  line-to-line  voltages  at  levels  that  approximate  the  experimental  results  in  Figure 
28.  Similarly,  part  (b)  shows  how  the  theoretical  input  voltages  are  transformed  through 
an  ideal  rectifier.  The  ideal  load  voltage  is  labeled  vi  ideal.  The  non-ideal  load  voltage  vi 
non-ideal  is  1.4V  lower  than  the  ideal  due  to  two  0.7V  losses  in  each  line-to-line  voltage 
incurred  by  the  diodes. 

The  experimental  results  of  the  physical  machine  show  approximately  a  three 
second  oscilloscope  capture  depicting  a  torque  impulse  that  accelerates  the  machine  and 
subsequently  allows  it  to  decay  to  nearly  zero  output.  These  results  can  be  seen  in  Figure 
27  and  Figure  28.  In  Figure  27  Ia  and  vl  fall  to  zero  while  Va  is  still  at  approximately  1.2 
volts  peak  {vpk).  This  is  expected  because  the  full  bridge  rectifier  suffers  from  voltage 
drop  across  two  diodes  in  each  phase.  The  diodes  in  this  circuit  were  measured  to  have 
roughly  0.6-0.7V  drop  each  in  the  forward  direction  which  accounts  for  the  missing 
1.2V. 
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(a) 


(b) 


Figure  26.  Theoretical  three-phase  diode  rectifier  operation  over  one  period,  (a) 
Sinusoidal  input  line  to  line  voltages,  (b)  Rectified  line  to  line 
voltages  and  load  voltage  with  ideal  and  non-ideal  diodes. 


Although  one  would  expect  to  see  the  rate  of  decay  decrease  after  the  diodes  stop 
conducting  and  the  machine  becomes  unloaded,  this  cannot  be  seen  in  the  results.  This  is 
likely  because  the  mechanical  damping  at  this  speed  is  much  larger  than  the  electrical 
loading  and  also  because  unloaded  generators  still  incur  some  losses. 
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Figure  27.  Results  of  a  single  input  torque  impulse  on  a  brushless  de  motor  with 

data  markers  indieating  the  eutoff  point  of  the  rectifier,  (a)  Voltage 
Vab  between  phases  a  and  c.  (b)  Current  ia  past  node  a.  (c)  Voltage  vl 

across  the  load. 


The  result  of  the  current  probe  being  used  between  phases  on  a  delta  connected 
machine  can  be  seen  in  Figure  28.  This  placement  manifests  itself  as  the  double  current 
spikes  seen  in  ia.  Also  visible  at  this  scale  is  the  six -pulse  ripple  output  of  vl  caused  by 
the  full  bridge  rectifier,  where  the  output  voltage  has  two  peaks  per  period  per  phase. 
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If  the  results  of  the  theoretical  six-pulse  rectifier  are  extended  and  only  the 
portions  that  were  experimentally  measured  are  plotted,  a  close  relationship  can  be  seen 
with  the  results  from  Figure  28.  This  correlation  can  be  seen  in  Figure  29. 


Figure  28.  Enlarged  selection  of  experimental  results  showing  double  current 
waveforms  and  six -pulse  rectified  voltage  from  Figure  27. 


Nomalized  Period 


Figure  29.  Theoretical  six -pulse  rectifier  results  representing  the  data  from 

Figure  28. 
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B. 


LINEAR  GENERATOR  SIMULINK  MODEL 


As  previously  stated,  the  model  generator  is  a  three-phase,  12-pole,  non- 
sinusoidal,  permanent-magnet,  linear  synchronous  machine.  The  model  parameters  are 
defined  for  motor  action  but  have  been  presented  in  the  result  plots  in  whichever  polarity 
is  most  intuitive.  The  workings  of  the  model  were  described  in  detail  in  Chapter  III  and 
the  full  MATLAB  and  Simulink  code  can  be  found  in  the  Appendix.  Several  adaptations 
of  the  model  have  been  created,  with  the  results  presented  here.  The  base  machine 
parameters  are  presented  in  Table  2.  These  are  the  machine  parameters  used  in  the  model 
in  each  scenario  unless  otherwise  stated. 

Plotting  the  portion  of  from  (17),  which  is  caused  by  the  permanent  magnets, 
allows  one  to  see  the  non-sinusoidal  distribution  of  the  machine,  and  is  shown  in  Figure 
30. 


Figure  30.  Magnetic  flux  profile  caused  by  permanent  magnets  as  experienced 

by  a  stator  winding  phase. 

The  exact  sizes  of  the  rotating  and  linear  machines  are  unknown  and  therefore 
some  assumptions  are  made  in  order  to  convert  the  model  from  the  original  design  of  a 
rotating  machine  into  the  linear  machine  that  is  being  simulating.  Specifically,  the  radius 
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of  the  rotor  in  the  rotating  maehine  and  the  length  of  the  travel  of  the  linear  maehine  are 
needed.  It  is,  therefore,  assumed  that  the  rotating  machine  has  a  rotor  radius  of  0.1 
meters.  Since  the  length  of  travel  of  a  linear  machine  must  be  equal  to  the  circumference 
of  the  rotating  machine,  the  travel  is  0.628  meters.  This  information  is  necessary  in  order 
to  convert  torque  to  force,  and  rotor  position  and  angular  velocity  to  translator  distance 
and  speed. 


Table  2.  Base  machine  parameters. 


4  =  24.5mH 

L  ^4.\2mH 

ms 

/I'  =0.15491F& 

m 

r  =9.1Q 

B  =0.0026Aw-5 

m 

J  -  0.0041kg 

a;  =-0.121 

K[  =  0.060 

a;  =  -0.009 

1,  Low  Frequency,  Low  Amplitude  Sinusoidal  Input  Force  Scenario 

In  the  sinusoidal  input  force  scenario  an  oscillating  input  is  being  considered, 
which  is  appropriate  to  linear  machines  associated  with  wave  energy  conversion.  While  a 
purely  sinusoidal  input  is  unrealistic  in  the  best  of  ocean  environments,  it  provides  a 
study  of  the  ideal  wave  energy  converter  input.  As  will  be  seen,  the  effects  of  the 
centering  spring  cause  a  non-sinusoidal  sum  of  forces  on  the  machine  that  will 
nevertheless  result  in  significant  harmonic  distortion. 

For  the  input  torque  to  the  model  generator,  a  sine  wave  with  an  amplitude  of 
0.8N-m  and  frequency  of  IHz  is  used.  Converting  torque  to  force,  where  F  -tR,  gives 
an  equivalent  sinusoidal  input  force  with  an  amplitude  of  8N,  still  at  the  same  frequency 
of  IHz.  This  is  seen  in  Figure  31,  where  the  torque  that  is  used  in  the  model  is  seen  in 
part  (a)  and  the  equivalent  force  is  seen  in  part  (b).  Taking  into  account  the  force  of  the 
centering  spring,  the  resulting  applied  mechanical  force  to  the  stator  is  given  in  Figure  32. 

This  input  causes  the  translator  to  oscillate  as  depicted  in  Figure  33,  where  part 
(a)  represents  the  velocity  (V)  of  the  oscillator  in  the  plane  of  movement  and  part  (b) 
depicts  the  position  of  the  translator  in  time. 
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While  the  mechanical  motion  of  the  translator  is  discussed  in  terms  of  velocity 
and  force,  for  electrical  variables  they  are  still  referred  to  by  the  angle-based  values  of 
0^^  and  0)^^.  The  response  of  these  variables  is  displayed  in  Figure  34,  where  can  be 
seen  in  part  (a)  in  rad/s.  The  translator  electrical  angle  9^^  covers  an  electrical  cycle 
(modulated  to  between  0-27i)  approximately  six  times  for  each  time  the  translator 
mechanically  travels  one  length  of  the  machine.  This  is  because  there  are  12  poles  (six 
pole-pairs)  in  the  machine. 


(a) 


(b) 


Figure  3 1 .  Sinusoidal  input  functions  over  three  periods,  (a)  Input  torque,  (b) 

Equivalent  input  force. 
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Figure  32.  Sinusoidal  input  force  and  resulting  applied  mechanical  force  to  the 
translator  taking  into  account  the  force  of  the  spring. 


(a) 


(b) 


Figure  33.  Translator  motion  over  three  periods,  (a)  Translator  velocity,  (b) 

Translator  position. 
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(a) 


Time  (s) 

Figure  34.  Electrical  motion,  (a)  Angular  velocity  .  (b)  Electrical  position  0^^ . 


The  stator  flux  linkage  from  (17),  is  depicted  in  Eigure  35,  where  the  system 

response  is  seen  in  (a)  and  a  more  detailed  view  of  the  waveforms  is  seen  in  (b).  Note  that 
the  waveforms  in  (b)  are  very  similar  to  those  in  Eigure  30,  with  the  differences  coming 
from  the  flux  produced  by  currents  in  the  stator  windings. 

From  (7)  it  can  be  seen  that  the  voltage  developed  in  the  stator  windings  is  in 
large  part  representative  of  the  rate  of  change  of  the  flux  linkage.  Comparing  the  shape  of 
the  waveforms  of  and  ,  this  relationship  can  be  seen  to  play  out,  where  the 
roughly  triangular  k^^^^  waveforms  produce  roughly  trapezoidal  waveforms.  The 
waveforms  for  can  be  seen  in  Figure  36  and  Figure  37.  The  stator  currents  are 

also  on  display  in  these  figures.  The  current  waveform  shape  is  caused  from  a 
combination  of  factors,  but  they  are  largely  controlled  by  the  effects  of  the  diode 
rectifiers  which  are  feeding  a  common  dc  bus. 
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Figure  35.  Stator  flux  linkages  with  the  reetangle  in  (a)  indicating  the 

selection  displayed  in  (b)  in  greater  detail. 


Time  (s) 


Figure  36.  Stator  machine  variables  over  three  periods,  (a)  Stator  voltages  . 

(b)  Stator  currents  .  Black  boxes  represent  data  selected  for 
display  in  higher  detail  in  Figure  37. 
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(a) 


Figure  37.  Stator  machine  variables  over  a  fraction  of  a  period  from  Figure  36. 

(a)  Stator  voltages  .  (b)  Stator  currents  . 


The  stator  voltages  are  rectified  and  applied  to  a  common  dc  bus,  which  can 

be  seen  in  Figure  38  and  Figure  39.  Only  the  bus  voltage  is  discussed  because  the  load  is 
a  purely  resistive  150Q  resistor.  The  dc  bus  has  a  single  l.lmF  capacitor  in  parallel  with 
the  load  resistor  in  order  to  help  smooth  the  output  voltage  to  this  load.  The  resultant  RC 
time  constant  (tkc),  given  by  -r-C,  equals  0.183s.  This  time  constant  affects  the  rate 
of  change  of  the  voltage  on  the  dc  bus  by 


v(t)  -  VqC  when  discharging, 


(32) 


and 


m- 


-t/  \ 


\-e' 


when  charging. 


(33) 


where  v(t)  is  the  voltage  at  time  t  and  vq  is  some  initial  voltage  at  ^  -  ^ .  The  practical 
effects  of  this  time  constant  can  be  seen  by  comparing  parts  (a)  and  (b)  in  Figure  39.  In 
the  first  portion  of  part  (a),  one  can  see  the  six -pulse  related  ripple  that  is  caused  by  the 
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individual  phases  seen  in  part  (b).  Note  that  even  this  ripple  is  smoother  than  the  spiky 
waveforms  of  part  (b).  While  this  is  an  effect  of  the  filter  capacitor  time  constant,  the 
more  dramatic  effect  is  seen  by  focusing  on  the  latter  portion  of  part  (b),  beginning 
around  1.5  seconds.  At  this  point,  the  voltage  ripple  fades  and  the  bus  voltage  is  riding  on 
the  capacitor.  Another  effect  of  the  filter  capacitor  is  the  apparent  phase  shift  between 
part  (a)  and  part  (b)  caused  by  the  capacitor  as  it  delays  the  voltage  rise  and  fall. 


Figure  38.  Voltage  on  the  dc  bus  over  three  periods  of  translator  motion. 

Power  of  the  machine  can  be  quantified  by  the  mechanical  power  in,  where 

P.  =  =  FV,  (34) 

or  by  the  electrical  power  out,  as 

P  =  vi,  (35) 

where  each  of  these  equations  provides  equivalent  units.  Using  (34)  and  (35)  the 
mechanical  power  put  into  the  machine  and  the  electrical  power  provided  to  the  load  can 
be  calculated. 
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Time(s) 

Figure  39.  Demonstration  of  the  filter  eapacitor’s  effeets.  (a)  Voltage  on  the  de 
bus  over  a  portion  of  a  period,  (b)  Currents  from  individual  stator 
phases  as  they  pass  onto  the  dc  bus  over  the  same  elapsed  time. 

To  find  the  total  energy  into  and  out  of  the  machine,  the  respective  power  over 
time  is  integrated.  By  comparing  these  two  values,  the  machine  efficiency  can  be 
uncovered.  These  calculations  have  been  performed  in  MATLAB  and  can  be  viewed  in 
the  Appendix.  It  should  be  noted  that  while  negative  power  is  defined  if  the  machine  is 
moving  one  direction  while  the  force  is  applied  in  the  other,  this  is  still  power  that  is 
being  applied  and  so  the  absolute  value  of  this  power  has  been  used  to  calculate  the 
machine  input  power  and  energy. 

The  results  for  this  scenario  are  that  mechanical  energy  into  the  machine  was 
18.46W-S  and  the  electrical  energy  delivered  to  the  load  was  9.98W-S  for  a  total 
efficiency  of  54.1%. 
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2,  Low  Frequency,  Low  Amplitude  Square  Input  Force  Scenario 

In  this  scenario  the  input  forcing  function  will  be  a  square  wave  at  IHz  with  an 
amplitude  of  ±0.565  N-m  of  torque  or  5.65N  of  foree.  This  value  is  the  rms  value  of  the 
amplitude  of  sine  input  force  in  order  to  make  them  more  comparable.  This  is  a  largely 
aeademie  scenario  since  a  square  input  funetion  at  this  frequency  is  not  realistic  to  any 
real  world  energy  sources  for  which  a  linear  generator  might  be  applied.  It  can  be  a  erude 
approximation  of  what  would  oeeur  in  a  free-piston  Stirling  engine  powered  generator  or 
perhaps  even  less  sophisticated,  in  a  linear  generator  masquerading  as  an  internal 
combustion  engine  cylinder,  although  both  of  these  applications  would  achieve  a  higher 
frequency.  These  higher  frequency  applications  will  be  the  topic  of  another  scenario. 

The  described  input  can  be  seen  in  Figure  40.  The  force  resulting  from  this  input 
foree  combined  with  the  force  of  the  spring  can  be  seen  in  Figure  41.  The  combined 
mechanical  force  acting  on  the  translator  is  more  intuitive  when  viewed  over  the  square 
input  function,  where  the  force  of  the  spring  adds  to  the  input  force  while  returning  to 
center  and  subtracts  as  the  machine  approaches  the  extremes  of  its  travel.  The  motion  of 
the  translator  is  depicted  in  Figure  42. 


(a) 


Figure  40.  Square  wave  input  functions  over  three  periods,  (a)  Input  torque,  (b) 

Equivalent  input  force. 
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Figure  41 .  Square  input  force  and  resulting  applied  mechanical  force  to  the 
translator  taking  into  account  the  force  of  the  spring. 


Figure  42.  Mechanical  motion  of  the  translator  resulting  from  a  square  input 

function. 
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The  flux  linkage  is  not  differentiable  from  the  previous  seenario  and  does  not 
merit  repeating.  The  stator  voltages  and  currents  are  very  similar  to  the  sine  input 
scenario,  except  that  the  harmonics  present  in  the  flux  linkage  are  more  exaggerated. 
These  are  seen  in  Figure  43  and  Figure  44,  with  Figure  44  part  (a)  clearly  showing  the 
harmonic  presence  where  the  voltage  wave  forms  deviate  from  trapezoidal. 


(a) 


0  0.5  1  1.5  2  2.5  3 

Time  (s) 

(b) 


Time  (s) 

Figure  43.  Stator  machine  variables  over  three  periods  of  the  square  input 
function,  (a)  Stator  voltages  .  (b)  Stator  currents  .  Black 
boxes  show  data  selected  for  display  in  greater  detail  in  Figure  44. 


The  dc  bus  voltages  behave  in  a  similar  manner  to  the  previous  scenario,  except 
that  the  more  abrupt  changes  in  speed  result  in  a  more  distinct  change  from  small  ripples 
in  the  bus  voltage  to  clearly  riding  on  the  capacitor  as  it  discharges  through  the  load 
according  to  its  time  constant.  This  is  portrayed  in  Figure  45  and  Figure  46. 
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(a) 


Figure  44.  Stator  machine  variables  over  a  fraction  of  a  period  from  Figure  43. 

(a)  Stator  voltages  .  (b)  Stator  currents  . 


O) 


Figure  45.  Voltage  on  the  dc  bus  over  three  periods  of  the  square  input  function. 
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Figure  46.  Demonstration  of  the  filter  capacitor’s  effects,  (a)  Voltage  on  the  dc 
bus  over  a  portion  of  a  period,  (b)  Currents  from  individual  stator 
phases  as  they  pass  onto  the  dc  bus  over  the  same  elapsed  time. 

To  determine  efficiency,  the  energy  into  and  out  of  the  machine  is  again 
calculated  as  in  the  previous  scenario.  Here  the  results  are  that  the  mechanical  energy  into 
the  machine  is  17.45W-S  while  the  energy  delivered  to  the  load  is  9.08W-S.  This  leads  to 
a  similar  efficiency  as  before  of  52.0%. 

3,  Floating  Frequency,  High  Amplitude  Impulse  Input  Force  Scenario 

It  is  apparent  through  experimenting  with  various  input  functions,  that  if  a  fixed 
frequency  alternating  impulse  is  applied  to  the  generator  the  machine  will  always 
stabilize  so  that  the  impulse  is  firing  while  the  translator  is  at  its  maximum  travel  and  in 
the  direction  to  extend  that  travel.  By  doing  so,  the  machine  would  always  find  its  least 
efficient  operating  point,  achieving  efficiencies  in  the  single  digits  at  best. 

In  order  to  simulate  a  crankless  internal  combustion  engine,  it  is  necessary  to  fire 
the  impulse  at  a  more  advantageous  point  in  the  translator  travel.  The  impulse  is  triggered 
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to  act  in  the  direction  of  motion  of  the  translator  as  the  translator  crosses  the  center  of  its 
travel  which  corresponds  to  zero  force  applied  by  the  spring.  In  this  way  the  maehine  will 
naturally  achieve  resonance  after  some  time.  This  is  quite  different  from  how  the  forces 
act  on  a  piston  in  an  internal  combustion  engine,  but  the  desired  takeaway  is  to 
understand  the  operation  of  a  linear  generator  that  is  not  frequency  controlled  and  that 
that  has  an  intermittent  input.  In  this  way,  this  scenario  explores  different  behaviors  in 
linear  generators  than  previous  results  and  is  still  in  the  vein  of  a  erankless  internal 
combustion  engine  driven  machine. 

For  the  results  shown,  the  input  impulse  force  is  applied  for  a  duration  of  one 
simulation  step,  which  is  20ps.  The  impulse  is  initially  ±50kN  and  quickly  ramps  down 
to  its  final  value  of  ±9kN.  The  high  initial  value  is  to  aid  the  machine  in  more  quickly 
reaching  a  steady  state  value.  The  steady  state  values  of  the  input  impulses  ean  be  seen  in 
Figure  47.  The  frequency  is  controlled  by  the  dynamics  of  the  machine;  it  settles  out  at 
roughly  70.4Hz.  If  you  equate  this  to  the  speed  of  a  piston  in  an  internal  eombustion 
engine,  the  engine  would  be  operating  at  4225RPM  which  falls  within  the  range  of  nearly 
all  automobile  engines.  Note  that  the  dimensions  that  have  been  assigned  to  this  machine 
are  not  compatible  with  the  dimensions  of  an  automobile  engine,  so  the  parallels  between 
them  should  be  filtered  accordingly. 


Figure  47.  Applied  input  impulses,  (a)  Impulse  as  torque,  (b)  Impulse  as  force. 
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In  Figure  48  the  external  forces  acting  on  the  translator,  including  the  input 
impulse  and  the  force  of  the  spring,  are  contrasted  with  the  motion  of  the  translator.  It  is 
apparent  that  the  translator  is  in  the  center  of  its  travel  when  the  impulse  is  applied.  This 
can  be  seen  by  the  position  of  the  translator  and  by  the  force  of  the  spring  passing 
through  zero.  Some  quantities  in  Figure  48  have  been  greatly  down-scaled  in  order  to 
more  intuitively  give  a  sense  of  the  operation  of  the  machine;  particular  attention  should 
be  paid  to  their  multipliers  in  the  legend. 


Figure  48.  External  forces  and  their  relationship  to  translator  motion. 

The  stator  voltages  and  currents  for  the  three  phases  can  be  seen  in  Figure  49.  The 
voltages  remain  trapezoidal  as  before;  the  ringing  that  is  visible  is  due  to  instability  in  the 
model  and  can  be  ignored.  The  current  waveforms  remain  largely  the  same  as  in  previous 
scenarios. 
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(a) 


(b) 


Figure  49.  Stator  variables  over  half  of  a  translator  eycle.  (a)  The  three  phases  of 
stator  voltage  .  (b)  The  three  phases  of  stator  current  . 

The  largest  change  is  seen  in  the  dc  bus  quality.  There  has  been  a  substantial 
increase  in  the  quality  of  the  power  on  the  dc  bus.  This  is  expected  at  higher  frequencies, 
where  the  ripple  can  be  more  easily  compensated  for  by  the  capacitor.  The  difference  in 
dc  bus  voltage  quality  can  been  seen  in  Figures  50  and  51.  The  average  dc  bus  voltage 
can  be  seen  to  build  up  and  then  reach  steady  state  at  about  50V.  This  is  much  different 
from  the  previous  scenarios  which  saw  the  dc  bus  voltage  in  a  periodic  state  of  flux  with 
large  transients  that  remained  during  steady  state.  The  ripple  that  is  visible  in  Figure  5 1  is 
no  longer  the  ripple  associated  with  the  six-pulse  rectifier  which  occur  six  times  each 
electrical  cycle,  but  rather  the  ripple  of  the  translator  moving  over  its  mechanical  cycle. 
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Figure  50.  Transient  into  steady  state  de  bus  voltage. 


Figure  5 1 .  Steady  state  de  bus  voltage  ripple. 


Here  it  ean  be  measured  that  the  total  de  bus  voltage  ripple  represents  only  2.7% 
of  the  de  bus  voltage.  In  the  sinusoidal  scenario  it  was  greater  than  50%  of  the  de  bus 
voltage  and  in  the  square  wave  scenario  it  was  roughly  44%. 
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The  efficiency  of  the  machine  in  the  scenario  was  calculated  after  the  machine 
had  reached  steady  state  to  avoid  the  acceleration  during  transient  portion.  The  final 
efficiency  number  was  71.2%.  When  the  input  impulse  was  administered  at  a  fixed 
frequency  the  observed  result  was  between  0.2-3%  efficiency. 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


Chapter  V  will  discuss  the  things  that  have  been  learned  from  this  process.  It  will 
also  discuss  the  shortcomings  of  the  model,  its  assumptions  and  its  omissions  of  actual 
phenomena.  In  discussing  these  faults,  ideas  for  future  work  to  address  them  will  be 
injected.  There  is  also  potential  for  future  work  in  imagining  new  ways  to  apply  the 
model  in  its  present  form. 

It  was  shown  that  a  model  for  a  three-phase  non-sinusoidal  permanent  magnet 
synchronous  machine  can  be  developed  in  Simulink  and  that  this  machine  can  be 
manipulated  to  represent  a  linear  machine.  Both  by  use  of  the  model  and  the  physical 
machine,  it  was  demonstrated  that  motors  and  generators  can  be  used  in  the  reverse  of  the 
operation  for  which  they  were  designed. 

Power  quality  and  efficiency  present  large  obstacles  to  linear  machines,  especially 
if  they  are  going  to  operate  at  low  frequencies.  There  is  a  need  to  explore  the  feasibility 
of  better  and  more  elaborate  output  filters  to  clean  up  the  power  quality  on  the  dc  bus. 
Where  efficiency  is  concerned,  there  is  still  a  great  deal  to  learn  and  explore,  even  within 
the  model  here.  There  are  dozens  of  variables  that  can  be  played  with  to  better  understand 
their  effects  on  efficiency,  including  the  spring  constant,  the  stator  resistance  ,  the 

leakage  inductance  ,  the  loss  coefficient  or  the  inertia  of  the  translator  J  . 

Major  challenges  in  working  with  these  variables  are  the  second  and  third  order 
effects.  A  change  to  one  variable  affects  the  overall  operation  of  the  machine  making 
comparisons  difficult.  For  instance,  a  change  in  J  could  make  the  machine  run  faster  and 
the  translator  travel  farther.  This  could  produce  more  power  but  at  lower  efficiency; 
increased  translator  travel  would  also  decrease  power  density  in  the  machine. 

To  be  truly  useful,  the  model  needs  to  take  into  account  both  cogging  torque 

(especially  at  low  speeds)  and  saturation.  This  goes  back  to  the  previous  point,  in  order  to 

quantify  any  results  the  state  of  the  machine  must  be  better  understood.  As  a  general  rule 

of  thumb,  electrical  machines  are  more  efficient  as  they  become  larger  and  faster.  The 

goodness  factor  says  this,  although  it  doesn’t  strictly  apply  to  permanent  magnet 
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machines.  This  is  problematic  if  efficiency  is  the  goal,  but  the  saturation  limit  of  the  core 
is  unknown. 
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APPENDIX  A.  GENERIC  SIMULINK  MODEL 


1,  Top  Level  View 
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2,  Permanent  Magnet  Machine 

a.  Top  Level  View 
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c.  Electrical  Torque 


i_Os  Product 
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d. 


Electrical  Power 


p_Os 


71 


3 


Ks  Transform 
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4 


Inverse  Ks  Transform 


5,  Balance  of  Forces  Equation 


To  Workspace4 
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6 


Diode  Rectifier 


a.  Top  Level  View 
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b. 


A-phase  Rectifier 
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APPENDIX  B.  LINEAR  GENERATOR  INITIAL  CONDITIONS  FILE 


clear 
close  all 
clc 


1,  PM  Machine  IC  Data 


poles  =  12;  %number  of  poles 
tstop=7 ; 
tstep=. 00002 ; 
rs=9.1; 

Lms=. 00412; 

Bm  =  0.0026; 

Ll s=24 . 5e-3 ; 

Lss=3/2*Lms+Ll s ; 

Inertia  =  .0041; 

Vdc=160 ; 

radius  =  0.1;  %  (m) 

fund_amp  =  0.6; 
fund_freq  =  2*pi ; 
harmonicl  =  5; 
harml_amp  =  .1; 
harmonic2  =  3; 
harm2_amp  =  .1; 

lambda_m  =  .1549; 

K3=-.121; 

K5=.060; 

K7=-.009; 


2,  Diode  Rectifier  IC  Data 


Ls=le-3;  %original 
Lsl=Ls/2 ; 

Ls2=Ls/2; 

Rs=le-3 ; 

Rload=50; 

Cdc=1100e-6; 

Published  with  MATLAB®  R2014a 
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APPENDIX  C.  LOW  FREQUENCY,  LOW  AMPLITUDE 
SINUSOIDAL  INPUT  FORCE  SCENARIO 


1,  Simulink  Input  Force 


2,  Defining  Variables 


%  Linear  Generator  Post  Run  File 
format  compact 
close  all  ; 
clc 

fig  =  0; 
ias=iabcs(:  ,1)  ; 
ibs=iabcs(: ,2) ; 
ics=iabcs(: ,3) ; 

vas=vabcs( : , 1) ; 
vbs=vabcs( : , 2) ; 
vcs=vabcs( : , 3) ; 

power  =  vdc.*idc; 

energy_out  =  sumCpower)/length(timev)*max(timev) 

energy_i n  =  sum(abs (i nput_torque . *omega_mech) ) /I ength(ti mev)*max(ti mev) 

efficiency  =  (energy_out  /  energy_in)  *  100 

funtime  =  linspaceC0,2,1000) ; 
fund  =  sin(2*pi*funtime) ; 
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harms  =  K3.*sin(2*pi*funtime*3)  ; 
harms  =  K5.*sin(2*pi*funtime*5)  ; 
harm?  =  K7.*sin(2*pi*funtime*7)  ; 

zNormal  =  theta_mech  ./  (2*pi);  %  normalized  stator  position 

z  =  ZNormal  .*  2*pi  *  radius;  %  translator  position  (m) 

V  =  omega_mech  .*  radius;  %  (m/s) 

inputForce  =  input_torque  /  radius;  %  (N) 
springForce  =  tau_spring  /  radius;  %  (N) 
electricForce  =  Te  /  radius;  %  (N) 

anerOB.  ■OUf 


BnertaFinM 


3,  Lambdas 


fig  =  fig  +1; 

figure  (fig)  ; 
subplot(2,l,2) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev,lambda_abc(: ,3) , . . . 
‘Linewidth’ ,2) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
title(‘(b)’) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

%  xlim([l,1.5]) 
xlim([1.3  1.4]) 

% 

subplot(2,l,l) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 

‘  Li  newi dth ’ , 1) 
title(‘  (a)  ’) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 

%  xlim([l,1.5]) 
xlim([0  3]) 

annotation(‘ rectangle’ , . . . 

[0.466625  0.585461689587424  0.024625  0 . 335952848722982] ,’ Li neWi dth ’, 2 ,.. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 
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(b) 


4,  Forces  at  play 


fig  =  fig+1; 
figure  (fig) 

plot(timev,-v,timev,inputForce+springForce, ’ Li newi  dth ’ ,2) 
hold  on 

plot(timev,-electricForce,  ’r’) 
hold  off 

legend(‘Velocity  (m/s)  ’  ,  ’  F_i_n+F_{spri ng}  (N)  ’  , ’F_{electric}  (N)  ’  ,  .  .  . 

‘ Locati on ’ , ’ Southeast ’ ) 
xlim([0  3]) 
ylabel  (‘Amplitude’) 

fig  =  fig+1; 
figure  (fig) 

pi ot (timev, input Force, timev,inputForce+sp ring Force, ’Li newi dth’ ,2) 
legend( ‘ F_i_n ’ , ’ F_i_n+F_{spri ng} ’ , ’ Location ’ , ’ Southwest ’ ) 
xlim([0  3]) 

xlabel(‘Time  (s)’);  yl  abel  ( ‘  Force  (N)’) 
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Force  (N)  Amplitude 
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5. 


Position  &  Velocity  (Mechanical) 


fig  =  fig+1; 
figure  (fig) 

plot(timev,z,timev,v, 'Linewidth' ,2) 
xlim([0  3]) 
ylabel (‘Amplitude’) 
xlabel(‘Time  (s)’) 

legend(‘Position  (m) ’ , 'Velocity  (m/s) ’ , ’Location’ , ’Northwest’) 


6,  ABC  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs)  ;  title(‘  (a)  ’) 

1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 
xlim([0  3]);  ylim([-12,12]) 
ylabel  (‘Voltage  (V)’);  xlabel(‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.460625  0.587426326129664  0.026875  0.332023575638502] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 
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subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics)  ;  title(‘  (b)  ’) 
legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
xlim([0  3]) 

ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.4625  0.113948919449902  0.025  0.332023575638507] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 


Time  (s) 


7,  ABC  Variables  short  timeframe 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs, ’Linewidth’ ,2) ;  title(‘ (a) ’) 
1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 
ylim([-12,12]) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics,  ’Linewidth’  ,2)  ;  title(‘  (b)  ’) 
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legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 


(a) 


(b) 


8,  qdO  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(3,l,l) 
plot(timev,i_qdos(:  ,1)) 
xlim([0  3]) 
legend(‘i_q’) 
subplot(3,l,2) 
plot(timev,i_qdos(:  ,2)) 
xlim([0  3]) 
legend(‘i_d’) 
subplot(3,l,3) 
plot(timev,i_qdos(:  ,3)) 
xlim([0  3]) 
legend(‘i_0’) 


iqs=i_qdos(: ,1) ; 
ids=i_qdos(: ,2) ; 
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iOs=i_qdos(: ,3) ; 
theta_r=i_qdos(: ,4) ; 
iqd_angle=angle(iqs-li*ids)*180/pi  ; 


9,  Velocity  and  Position  (Electrical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,wre_data, 'Linewidth' ,2) 
titleC'  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Angular  Velocity  (rad/s)’) 
xlim([0  3]);  ylim([-125  125]) 

%1 egend ( ‘ \omega_r_e ’ ) 

subplot(2,l,2) 

plot(timev, theta_elec, ’ r ’ , ' Li newi dth ’ , 2) 
title(‘(b)’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Electrical  Position  (rad)’) 
xlim([0  3]) 

%1 egend ( ‘ \theta_r_e ’ ) 
ylim([0,2*pi]) 
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(a) 


Time  (s) 


10,  Velocity  and  Position  (Mechanical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 
plot(timev,  V,  ’  Li  newi  dth  ’  ,  2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Velocity  (m/s)’) 
xlim([0  3]); 

%1 egend ( ‘ Vel oci ty ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

pi  ot (ti  mev , z , ’ r ’ , ’ Li newi dth  ’  ,2) 
title(‘(b)’) 

xlim([0  3]);  ylim([0  .7]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Position  (m)’) 

%1  egend (‘Position’) 
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(a) 


11,  DC  Variables 


fig  =  fig+1; 
figure  (fig) 

pi ot (ti mev , vdc , ’Li newi dth ’ ,2) 

xlim([0  3]);  ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,  vdc,  ’  Li  newi  dth  ’  ,  2)  ;  ti  tl  e(  ‘  (a)  ’) 

xlim([1.25  1.65]):  %ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

subplot(2,l,2) 
figure  (fig) 

plot(timev,idcabc(:  ,1)  ,timev,idcabc(:  ,2)  ,timev,idcabc(:  ,  3))  ; 
title(‘(b)’) 

legend(‘i_a  dc’,’i_b  dc’,’i_c  dc’ , ’Location’ , ’Best’) 
xlim([1.25  1.65]):  %ylim([-l  17]) 
xlabel  (‘Time(s) ’)  ;  ylabel  (‘Current  (A)’) 
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Current  (A)  Voltage  (V)  Voltage  (V) 


Time  (s) 

(a) 


(b) 


Time(s) 
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12.  Lambda  Harmonics 


fig  =  fig+1; 
figure  (fig) 

plot(funtime,0.1549*(fund+harm3+harm5+harm7) , ’Linewidth’ ,2) 
xlim([0  max(funtime)] ) 

%1 egend ( ‘ \1 ambda_P_M ’ ) 

xlabel  (‘Time  (arbitrary)’);  ylabel  (‘Magnetic  Flux  (wb)’) 


13,  Input  Torque  and  Force 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,input_torque, ’ Linewidth’ ,2) 

title(‘(a)’) 

xlim([0  3]); 

%legend(‘\tau_i_n’ , ’Location’ , ’Southeast’) 
ylabel (‘Torque  (N\cdotm)’) 
xlabel  (‘Time  (s)’) 
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subplot(2,l,2) 

plot(timev, i nput Force, ’ r’ , ' Li newi dth ’ , 2) 

titleCCb)’) 

xlim([0  3]) 

%legendC‘F_i_n’ , ’Location' , ’Southeast’) 
ylabel  (‘ Force  (N)’) 
xlabel  (‘Time  (s)’) 


(a) 


(b) 


Published  with  MATLAB®  R2014a 
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APPENDIX  D.  LOW  FREQUENCY,  LOW  AMPLITUDE  SINE 
INPUT  WITH  MULTIPLE  HARMONICS  ADDED 


1,  Simulink  Input  Force 


2,  Defining  Variables 


%  Linear  Generator  Post  Run  File 
format  compact 
close  all  ; 
clc 

fig  =  0; 
ias=iabcs(:  ,1)  ; 
ibs=iabcs(: ,2) ; 
ics=iabcs(: ,3) ; 

vas=vabcs( : , 1) ; 
vbs=vabcs( : , 2) ; 
vcs=vabcs( : , 3) ; 

power  =  vdc.*idc; 

energy_out  =  sumCpower)/length(timev)*max(timev) 

energy_i n  =  sum(abs (i nput_torque . *omega_mech) ) /I ength(ti mev)*max(ti mev) 

efficiency  =  (energy_out  /  energy_in)  *  100 

funtime  =  linspaceC0,2,1000) ; 
fund  =  sin(2*pi*funtime) ; 
harm3  =  K3.*sinC2*pi*funtime*3) ; 
harms  =  K5.*sinC2*pi*funtime*5) ; 
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harm?  =  K7.*sin(2*pi*funtime*7)  ; 


zNormal  =  theta_mech  ./  (2*pi);  %  normalized  stator  position 

z  =  ZNormal  .*  2*pi  *  radius;  %  translator  position  (m) 

V  =  omega_mech  .*  radius;  %  (m/s) 

inputForce  =  input_torque  /  radius;  %  (N) 
springForce  =  tau_spring  /  radius;  %  (N) 
electricForce  =  Te  /  radius;  %  (N) 


Bneratfcouf 


3,  Lambdas 


fig  =  fig  +1; 

figure  (fig)  ; 
subplot(2,l,2) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 

‘ Li newi dth ’ , 2) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
title(‘(b)’) 

xlabel(‘Time  (s)’);  yl abel  ( ‘Magnetic  Flux  (wb)’) 

%  xlim([l,l.  5]) 
xlim([1.3  1.4]) 

% 

subplot(2,l,l) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 

‘  Li  newi dth ’ , 1) 
title(‘  (a)  ’) 

xl  abel  (‘Time  (s)’);  yl  abel  (‘Magnetic  Flux  (wb)’) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 

%  xlim([l,l.  5]) 
xlim([0  3]) 

annotation(‘ rectangle’ , . . . 

[0.466625  0.585461689587424  0.024625  0 . 335952848722982] ,’ Li neWi dth ’, 2 ,.. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 
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(a) 


(b) 


4,  Forces  at  play 


fig  =  fig+1; 
figure  (fig) 

plot(timev,-v,timev,inputForce+springForce, ’ Li newi  dth ’ ,2) 
hold  on 

plot(timev,-electricForce,  ’r’) 
hold  off 

legend(‘Velocity  (m/s)  ’  ,  ’  F_i_n+F_{spri ng}  (N)  ’  , ’F_{electric}  (N)  ’  ,  .  .  . 

‘ Locati on ’ , ’ Southeast ’ ) 
xlim([0  3]) 
ylabel  (‘Amplitude’) 

fig  =  fig+1; 
figure  (fig) 

pi ot (timev, input Force, timev,inputForce+sp ring Force, ’Li newi dth’ ,2) 
legend( ‘ F_i_n ’ , ’ F_i_n+F_{spri ng} ’ , ’ Location ’ , ’ Southwest ’ ) 
xlim([0  3]) 

xlabel(‘Time  (s)’);  yl  abel  ( ‘  Force  (N)’) 
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Force  (N)  Amplitude 


98 


5. 


Position  &  Velocity  (Mechanical) 


fig  =  fig+1; 
figure  (fig) 

plot(timev,z,timev,v, 'Linewidth' ,2) 
xlim([0  3]) 
ylabel (‘Amplitude’) 
xlabel(‘Time  (s)’) 

legend(‘Position  (m) ’ , 'Velocity  (m/s) ’ , ’Location’ , ’Northwest’) 


6,  ABC  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs)  ;  title(‘  (a)  ’) 

1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 
xlim([0  3]);  ylim([-12,12]) 
ylabel  (‘Voltage  (V)’);  xlabel(‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.460625  0.587426326129664  0.026875  0.332023575638502] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 
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subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics)  ;  title(‘  (b)  ’) 
legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
xlim([0  3]) 

ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.4625  0.113948919449902  0.025  0.332023575638507] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 


Time  (s) 

7,  ABC  Variables  short  timeframe 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs, ’Linewidth’ ,2) ;  title(‘ (a) ’) 
1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 
ylim([-12,12]) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics,  ’Linewidth’  ,2)  ;  title(‘  (b)  ’) 
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legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 


(a) 


Time  (s) 

(b) 


8,  qdO  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(3,l,l) 
plot(timev,i_qdos(:  ,1)) 
xlim([0  3]) 
legend(‘i_q’) 
subplot(3,l,2) 
plot(timev,i_qdos(:  ,2)) 
xlim([0  3]) 
legend(‘i_d’) 
subplot(3,l,3) 
plot(timev,i_qdos(:  ,3)) 
xlim([0  3]) 
legend(‘i_0’) 


iqs=i_qdos(: ,1) ; 
ids=i_qdos(: ,2) ; 
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iOs=i_qdos(: ,3) ; 
theta_r=i_qdos(: ,4) ; 
iqd_angle=angle(iqs-li*ids)*180/pi  ; 


1 

0 

-1 


"TTTT’^ — 


0.5 


1.5 


2.5 


9,  Velocity  and  Position  (Electrical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,wre_data, 'Linewidth' ,2) 
titleC'  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Angular  Velocity  (rad/s)’) 
xlim([0  3]);  ylim([-125  125]) 

%1 egend ( ‘ \omega_r_e ’ ) 

subplot(2,l,2) 

plot(timev, theta_elec, ’ r ’ , ' Li newi dth ’ , 2) 
title(‘(b)’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Electrical  Position  (rad)’) 
xlim([0  3]) 

%1 egend ( ‘ \theta_r_e ’ ) 
ylim([0,2*pi]) 
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(a) 


(b) 


10,  Velocity  and  Position  (Mechanical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 
plot(timev,  V,  ’  Li  newi  dth  ’  ,  2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Velocity  (m/s)’) 
xlim([0  3]); 

%1 egend ( ‘ Vel oci ty ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

pi  ot (ti  mev , z , ’ r ’ , ’ Li newi dth  ’  ,2) 
title(‘(b)’) 

xlim([0  3]);  ylim([0  .7]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Position  (m)’) 

%1  egend (‘Position’) 


103 


(a) 


11,  DC  Variables 


fig  =  fig+1; 
figure  (fig) 

pi ot (ti mev , vdc , ’Li newi dth ’ ,2) 

xlim([0  3]);  ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,  vdc,  ’  Li  newi  dth  ’  ,  2)  ;  ti  tl  e(  ‘  (a)  ’) 

xlim([1.25  1.65]):  %ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

subplot(2,l,2) 
figure  (fig) 

plot(timev,idcabc(:  ,1)  ,timev,idcabc(:  ,2)  ,timev,idcabc(:  ,  3))  ; 
title(‘(b)’) 

legend(‘i_a  dc’,’i_b  dc’,’i_c  dc’ , ’Location’ , ’Best’) 
xlim([1.25  1.65]):  %ylim([-l  17]) 
xlabel  (‘Time(s) ’)  ;  ylabel  (‘Current  (A)’) 


1U4 


Current  (A)  Voltage  (V)  Voltage  (V) 


Time  (s) 

(a) 


(b) 


Time(s) 
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12.  Lambda  Harmonics 


fig  =  fig+1; 
figure  (fig) 

plot(funtime,0.1549*(fund+harm3+harm5+harm7) , ’Linewidth’ ,2) 
xlim([0  max(funtime)] ) 

%1 egend ( ‘ \1 ambda_P_M ’ ) 

xlabel  (‘Time  (arbitrary)’);  ylabel  (‘Magnetic  Flux  (wb)’) 


13,  Input  Torque  and  Force 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,input_torque, ’ Linewidth’ ,2) 

title(‘(a)’) 

xlim([0  3]); 

%legend(‘\tau_i_n’ , ’Location’ , ’Southeast’) 
ylabel (‘Torque  (N\cdotm)’) 
xlabel  (‘Time  (s)’) 
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subplot(2,l,2) 

plot(timev, i nput Force, ’ r’ , ' Li newi dth ’ , 2) 

titleCCb)’) 

xlim([0  3]) 

%legendC‘F_i_n’ , ’Location' , ’Southeast’) 
ylabel  (‘ Force  (N)’) 
xlabel  (‘Time  (s)’) 


(a) 


(b) 
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APPENDIX  E.  LOW  FREQUENCY,  LOW  AMPLITUDE  SQUARE 

INPUT  FORCE  SCENARIO 


1,  Simulink  Input  Force 


2,  Defining  Variables 

%  Linear  Generator  Post  Run  File 
format  compact 
close  all  ; 
clc 

fig  =  0; 
ias=iabcs(:  ,1)  ; 
ibs=iabcs(: ,2) ; 
ics=iabcs(: ,3) ; 

vas=vabcs( : , 1) ; 
vbs=vabcs( : , 2) ; 
vcs=vabcs( : , 3) ; 

power  =  vdc.*idc; 

energy_out  =  sum(power)/length(timev)*max(timev) 

energy_i n  =  sum(abs (i nput_torque . *omega_mech) ) /I ength(ti mev)*max(ti mev) 
efficiency  =  (energy_out  /  energy_in)  *  100 


funtime  =  linspace(0,2,1000) ; 
fund  =  sin(2*pi*funtime) ; 
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harms  =  K3.*sin(2*pi*funtime*3)  ; 
harms  =  K5.*sin(2*pi*funtime*5)  ; 
harm?  =  K7.*sin(2*pi*funtime*7)  ; 

zNormal  =  theta_mech  ./  (2*pi);  %  normalized  stator  position 

z  =  ZNormal  .*  2*pi  *  radius;  %  translator  position  (m) 

V  =  omega_mech  .*  radius;  %  (m/s) 

inputForce  =  input_torque  /  radius;  %  (N) 
springForce  =  tau_spring  /  radius;  %  (N) 
electricForce  =  Te  /  radius;  %  (N) 


ttneratfcouf 


3,  Lambdas 


fig  =  fig  +1; 

figure  (fig)  ; 
subplot(2,l,2) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 
‘Linewidth’ ,2) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
title(‘(b)’) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

%  xlim([l,1.5]) 
xlim([1.3  1.4]) 

% 

subplot(2,l,l) 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 

‘  Li  newi dth ’ , 1) 
title(‘  (a)  ’) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 

%  xlim([l,1.5]) 
xlim([0  3]) 

annotation(‘ rectangle’ , . . . 

[0.466625  0.585461689587424  0.024625  0 . 335952848722982] ,’ Li neWi dth ’, 2 ,.. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 


no 


(a) 


(b) 


4,  Forces  at  play 


fig  =  fig+1; 
figure  (fig) 

plot(timev,-v,timev,inputForce+springForce, ’ Li newi  dth ’ ,2) 
hold  on 

plot(timev,-electricForce,  ’r’) 
hold  off 

legend(‘Velocity  (m/s)  ’  ,  ’  F_i_n+F_{spri ng}  (N)  ’  , ’F_{electric}  (N)  ’  ,  .  .  . 

‘ Locati on ’ , ’ Southeast ’ ) 
xlim([0  3]) 
ylabel  (‘Amplitude’) 

fig  =  fig+1; 
figure  (fig) 

pi ot (timev, input Force, timev,inputForce+sp ring Force, ’Li newi dth’ ,2) 
legend( ‘ F_i_n ’ , ’ F_i_n+F_{spri ng} ’ , ’ Location ’ , ’ Southwest ’ ) 
xlim([0  3]) 

xlabel(‘Time  (s)’);  yl  abel  ( ‘  Force  (N)’) 
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Force  (N)  Amplitude 
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5. 


Position  &  Velocity  (Mechanical) 


fig  =  fig+1; 
figure  (fig) 

plot(timev,z,timev,v, 'Linewidth' ,2) 
xlim([0  3]) 
ylabel (‘Amplitude’) 
xlabel(‘Time  (s)’) 

legend(‘Position  (m) ’ , 'Velocity  (m/s) ’ , ’Location’ , ’Northwest’) 


6,  ABC  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs)  ;  title(‘  (a)  ’) 

1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati  on ’ , ’ Southwest ’ ) 
xlim([0  3]);  ylim([-12,12]) 
ylabel  (‘Voltage  (V)’);  xlabel(‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.460625  0.587426326129664  0.026875  0.332023575638502] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 
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subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics)  ;  title(‘  (b)  ’) 
legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
xlim([0  3]) 

ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
annotation(‘ rectangle’ , . . . 

[0.4625  0.113948919449902  0.025  0.332023575638507] , ’Linewidth’ ,2, .. . 
‘ FaceCol or ’ , ’ f 1  at ’ ) ; 


Time  (s) 

7,  ABC  Variables  short  timeframe 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs, ’Linewidth’ ,2) ;  title(‘ (a) ’) 
1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 
ylim([-12,12]) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics,  ’Linewidth’  ,2)  ;  title(‘  (b)  ’) 


114 


legend (‘i_a_s’ , ’i_b_s’ , ’i_c_s' , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
xlim([1.25  1.35]) 


(a) 


(b) 


8,  qdO  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(3,l,l) 
plot(timev,i_qdos(:  ,1)) 
xlim([0  3]) 
legend(‘i_q’) 
subplot(3,l,2) 
plot(timev,i_qdos(:  ,2)) 
xlim([0  3]) 
legend(‘i_d’) 
subplot(3,l,3) 
plot(timev,i_qdos(:  ,3)) 
xlim([0  3]) 
legend(‘i_0’) 


iqs=i_qdos(: ,1) ; 
ids=i_qdos(: ,2) ; 
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iOs=i_qdos(: ,3) ; 
theta_r=i_qdos(: ,4) ; 
iqd_angle=angle(iqs-li*ids)*180/pi  ; 


9,  Velocity  and  Position  (Electrical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,wre_data, 'Linewidth' ,2) 
titleC'  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Angular  Velocity  (rad/s)’) 
xlim([0  3]);  ylim([-125  125]) 

%1 egend ( ‘ \omega_r_e ’ ) 

subplot(2,l,2) 

plot(timev, theta_elec, ’ r ’ , ' Li newi dth ’ , 2) 
title(‘(b)’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Electrical  Position  (rad)’) 
xlim([0  3]) 

%1 egend ( ‘ \theta_r_e ’ ) 
ylim([0,2*pi]) 
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(a) 


Time  (s) 


10,  Velocity  and  Position  (Mechanical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 
plot(timev,  V,  ’  Li  newi  dth  ’  ,  2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Velocity  (m/s)’) 
xlim([0  3]); 

%1 egend ( ‘ Vel oci ty ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

pi  ot (ti  mev , z , ’ r ’ , ’ Li newi dth  ’  ,2) 
title(‘(b)’) 

xlim([0  3]);  ylim([0  .7]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Position  (m)’) 

%1  egend (‘Position’) 
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(a) 


11,  DC  Variables 


fig  =  fig+1; 
figure  (fig) 

pi ot (ti mev , vdc , ’Li newi dth ’ ,2) 

xlim([0  3]);  ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,  vdc,  ’  Li  newi  dth  ’  ,  2)  ;  ti  tl  e(  ‘  (a)  ’) 

xlim([1.25  1.65]):  %ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

subplot(2,l,2) 
figure  (fig) 

plot(timev,idcabc(:  ,1)  ,timev,idcabc(:  ,2)  ,timev,idcabc(:  ,  3))  ; 
title(‘(b)’) 

legend(‘i_a  dc’,’i_b  dc’,’i_c  dc’ , ’Location’ , ’Best’) 
xlim([1.25  1.65]):  %ylim([-l  17]) 
xlabel  (‘Time(s) ’)  ;  ylabel  (‘Current  (A)’) 
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Time  (s) 


(a) 


(b) 
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12.  Lambda  Harmonics 


fig  =  fig+1; 
figure  (fig) 

plot(funtime,0.1549*(fund+harm3+harm5+harm7) , ’Linewidth’ ,2) 
xlim([0  max(funtime)] ) 

%1 egend ( ‘ \1 ambda_P_M ’ ) 

xlabel  (‘Time  (arbitrary)’);  ylabel  (‘Magnetic  Flux  (wb)’) 


13,  Input  Torque  and  Force 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,input_torque, ’ Linewidth’ ,2) 

title(‘(a)’) 

xlim([0  3]); 

%legend(‘\tau_i_n’ , ’Location’ , ’Southeast’) 
ylabel (‘Torque  (N\cdotm)’) 
xlabel  (‘Time  (s)’) 
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subplot(2,l,2) 

plot(timev, i nput Force, ’ r’ , ' Li newi dth ’ , 2) 

titleCCb)’) 

xlim([0  3]) 

%legendC‘F_i_n’ , ’Location' , ’Southeast’) 
ylabel  (‘ Force  (N)’) 
xlabel  (‘Time  (s)’) 


(a) 


(b) 
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APPENDIX  F.  HIGH  FREQUENCY,  HIGH  AMPLITUDE  SQUARE 

INPUT  FORCE  SCENARIO 


1,  Simulink  Input  Force 


2,  Defining  Variables 


format  compact 
close  all  ; 
clc 

fig  =  0; 
ias=iabcs(:  ,1)  ; 
ibs=iabcs(: ,2) ; 
ics=iabcs(: ,3) ; 

vas=vabcs( : , 1) ; 
vbs=vabcs( : , 2) ; 
vcs=vabcs( : , 3) ; 

power  =  vdc.*idc; 

energy_sum  =  (sum(abs(input_torque.*omega_mech))-sum(power) . . . 
-sum(Bm*omega_mech) ) /I ength (ti mev) 

energy_out  =  sum(power(length(ti mev) -50000 : length (ti mev))) /50000 
energy_in  =  sum(abs(input_torque(length(timev)-50000:length(timev)) . . . 
. *omega_mechCl ength (ti mev) -50000:1  ength (ti mev) )))/50000 
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efficiency  =  (energy_out  /  energy_in)  *  100 


funtime  =  1inspace(0,2,1000) ; 
fund  =  sin(2*pi*funtime) ; 
harm!  =  K3.*sin(2*pi*funtime*3)  ; 
harms  =  K5.*sin(2*pi*funtime*5)  ; 
harm?  =  K7.*sin(2*pi*funtime*7)  ; 

zNormal  =  theta_mech  ./  (2*pi);  %  normalized  stator  position 

z  =  ZNormal  .*  2*pi  *  radius;  %  translator  position  (m) 

V  =  omega_mech  .*  radius;  %  (m/s) 
tau_frict  =  omega_mech*.0026; 


frictionForce  =  tau_frict  /  radius; 
inputForce  =  input_torque  /  radius;  %  (N) 
springForce  =  tau_spring  /  radius;  %  (N) 
electricForce  =  Te  /  radius;  %  (N) 


ftnercW-  .sum  M 

BneroM^out^ 


e+IS 


3,  Lambda 


fig  =  fig  +1; 
figure  (fig)  ; 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 
‘  Li  newi dth ’ , 1) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 

%  xlim([l,1.5]) 
xlim([4.9  5]) 
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4,  Forces  at  play 


fig  =  fig+1; 
figure  (fig) 

plot(timev,-inputForce/10000,timev,-(inputForce+springForce)/10000, . . . 

timev,z,timev,v/100, ’ Linewidth’ ,2) 
legend (‘F_i_n/10  (kN) ’ , ’ F_i_n+F_{spri ng}/10  (kN) ’ , . . . 

‘Position  (m) ’Velocity/100  (m/s) ’Location’ , ’Southwest’) 

xl i m ( [max (ti mev) - . 1  max ( t i mev) ] ) 
xlabel(‘Time  (s)’);  ylabel  (‘Amplitude’) 
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5,  ABC  Variables 


fig  =  fig+1; 
figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs) ;  title(‘ (a) ’) 

1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics)  ;  title(‘  (b)  ’) 

1 egend (‘i_a_s’ , ’i_b_s’ , ’i_c_s’ , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
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Time  (s) 


6,  ABC  Variables  short  timeframe 


fig  =  fig+1; 
figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs, ’Linewidth’ ,2) ;  title(‘ (a) ’) 
1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 
xl i m ( [max (ti mev) - . 1  max ( t i mev) ] ) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics,  ’Linewidth’  ,2)  ;  title(‘  (b)  ’) 
1 egend (‘i_a_s’ , ’i_b_s’ , ’i_c_s’ , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
xlim([max(timev)-.l  max(timev)]) 
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Time  (s) 

(b) 

10  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


Time  (s) 


7,  qdO  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(3,l,l) 
plot(timev,i_qdos(:  ,1)) 
xlim([4.9  5]) 
legend(‘i_q’) 
subplot(3,l,2) 
plot(timev,i_qdosC:  ,2)) 
xlim([4.9  5]) 
legend(‘i_d’) 
subplot(3,l,3) 
plot(timev,i_qdos(:  ,3)) 
xlim([4.9  5]) 
legend(‘i_0’) 

iqs=i_qdos(: ,1) ; 

ids=i_qdos(: ,2) ; 

iOs=i_qdos( : , 3) ; 

theta_r=i_qdos ( : , 4) ; 

i qd_angl e=angl e (i qs-li *i ds) *180/pi ; 
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8,  Velocity  and  Position  (Electrical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,wre_data, ’Linewidth’ ,2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Angular  Velocity  (rad/s)’) 
xlim([4.9  5]);  %ylim([-125  125]) 

%1 egend ( ‘ \omega_r_e ’ ) 

subplot(2,l,2) 

plot(timev, theta_elec, ’ r ’ , ’ Li newi dth ’ , 2) 
title(‘(b)’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Electrical  Position  (rad)’) 
xlim([4.9  5]) 

%1 egend ( ‘ \theta_r_e ’ ) 

%ylim([0,2*pi]) 
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(a) 


Time  (s) 


9,  Velocity  and  Position  (Mechanical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 
plot(timev,  V,  ’  Li  newi  dth  ’  ,  2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Velocity  (m/s)’) 
xlim([4.9  5]); 

%1 egend ( ‘ Vel oci ty ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

pi  ot (ti  mev , z , ’ r ’ , ’ Li newi dth  ’  ,2) 
title(‘(b)’) 

xlim([4.9  5]);  ylim([0  2*pi/10]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Position  (m)’) 

%1  egend (‘Position’) 
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(a) 


10,  DC  Variables 


fig  =  fig+1; 
figure  (fig) 

pi ot (ti mev , vdc , ’Li newi dth ’ ,2) 

xlim([4.9  5]);  %ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fig  =  fig+1; 

figure  (fig) 

plot(timev,  vdc,  ’  Li  newi  dth  ’ ,  2)  ;  %title(  ‘  (a)  ’) 

%xlim([0.5  0.6]);  %%ylim([0  11]) 
xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fi  g=fi  g+1 
figure  (fig) 

plot(timev,idcabc(:  ,1)  ,timev,idcabc(:  ,2)  ,timev,idcabc(:  ,  3))  ; 
%title(‘(b)’) 

legend(‘i_a  dc’,’i_b  dc’,’i_c  dc’ , ’Location’ , ’Best’) 
xlim([2.34  2.367]);  %%ylim([-l  17]) 
xlabel  (‘Time(s) ’)  ;  ylabel  (‘Current  (A)’) 
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%  %%  Lambda  Harmonics 
%  fig  =  fig+1; 

% 

%  figure(fig) 

%  plot(funtime,0.1549*(fund+harm3+harm5+harm7) , ’Linewidth’ ,2) 
%  xlim([0  max(funtime)] ) 

%  %legend(‘\lambda_P_M’) 

%  xlabel(‘Time  (arbitrary)');  ylabel  (‘Magnetic  Flux  (wb)’) 

fig  = 

10 
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Time(s) 
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11,  Input  Torque  and  Force 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,input_torque, ’ Linewidth’ ,2) 

titleCCa)’) 

xlim([4.9  5]); 

%legend(‘\tau_i_n’ , ’Location’ , ’Southeast’) 
ylabel (‘Torque  (N\cdotm) ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

plot(timev, i nput Force, ’ r ’ , ’ Li newi dth ’ , 2) 

title(‘(b)’) 

xlim([4.9  5]) 

%legend(‘F_i_n’ , ’Location’ , ’Southeast’) 
ylabel  (‘ Force  (N)’) 
xlabel  (‘Time  (s)’) 


(a) 


Time  (s) 

(b) 


Time  (s) 
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APPENDIX  G.  FLOATING  FREQUENCY,  HIGH  AMPLITUDE 
IMPULSE  INPUT  FORCE  SCENARIO 


1,  Simulink  Input  Force 


U<  0 
&  NOT 
U/z  <  0 


Detect  Fall 
Negative 


U>  0 
&  NOT 
U/z>0 


Detect  Rise 
Positive 


n 


— ► 

Product! 


Step  Rale  Limiter 


2,  Defining  Variables 


format  compact 
close  all  ; 
clc 

fig  =  0; 
ias=iabcs(:  ,1)  ; 
ibs=iabcs(: ,2) ; 
ics=iabcs(: ,3) ; 

vas=vabcs( : , 1) ; 
vbs=vabcs( : , 2) ; 
vcs=vabcs( : , 3) ; 

power  =  vdc.*idc; 

energy_sum  =  (sum(abs(input_torque.*omega_mech))-sum(power) . . . 
-sum(Bm*omega_mech) ) /I ength (ti mev) 

energy_out  =  sum(power(length(ti mev) -50000 : length (ti mev))) /50000 
energy_in  =  sum(abs(input_torque(length(timev)-50000:length(timev)) . . . 
. *omega_mech(l ength (ti mev) -50000:1  ength (ti  mev) )))/50000 

efficiency  =  (energy_out  /  energy_in)  *  100 


-►  nput_torque 


To  Woricspace 
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funtime  =  1inspace(0,2,1000) ; 
fund  =  sin(2*pi*funtime) ; 
harms  =  K3.*sin(2*pi*funtime*3) 
harms  =  K5.*sin(2*pi*funtime*5) 
harm?  =  K7.*sin(2*pi*f  untime*?) 


zNormal  =  theta_mech  ./  (2*pi);  %  normalized  stator  position 

z  =  ZNormal  .*  2*pi  *  radius;  %  translator  position  (m) 

V  =  omega_mech  .*  radius;  %  (m/s) 
tau_frict  =  omega_mech*.0026; 

frictionForce  =  tau_frict  /  radius; 
inputForce  =  input_torque  /  radius;  %  (N) 
springForce  =  tau_spring  /  radius;  %  (N) 
electricForce  =  Te  /  radius;  %  (N) 


3,  Lambdas 


fig  =  fig  +1; 
figure  (fig)  ; 

plot(timev,lambda_abc(: ,1) ,timev,lambda_abc(: ,2) ,timev, lambda_abc(: ,3) , . . . 
‘  Li  newi dth ’ , 1) 

xlabel(‘Time  (s)’);  ylabel  (‘Magnetic  Flux  (wb)’) 

1 egend ( ‘ \1 ambda_a_s ’ , ’ \1 ambda_b_s ’ , ’ \1 ambda_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 

%  xlim([l,1.5]) 
xlim([4.9  5]) 
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4,  Forces  at  play 


fig  =  fig+1; 
figure  (fig) 

plot(timev,-inputForce/10000,timev,-(inputForce+springForce)/10000, . . . 

timev,z,timev,v/100, ’ Linewidth’ ,2) 
legend (‘F_i_n/10  (kN) ’ , ’ F_i_n+F_{spri ng}/10  (kN) ’ , . . . 

‘Position  (m) ’Velocity/100  (m/s) ’Location’ , ’Southwest’) 

xl i m ( [max (ti mev) - . 1  max ( t i mev) ] ) 
xlabel(‘Time  (s)’);  ylabel  (‘Amplitude’) 
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5,  ABC  Variables 


fig  =  fig+1; 
figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs) ;  title(‘ (a) ’) 

1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics)  ;  title(‘  (b)  ’) 

1 egend (‘i_a_s’ , ’i_b_s’ , ’i_c_s’ , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
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(a) 


Time  (s) 

6,  ABC  Variables  short  timeframe 


fig  =  fig+1; 
figure  (fig) 
subplot(2,l,l) 

plot(timev,vas,timev,vbs,timev,vcs, ’Linewidth’ ,2) ;  title(‘ (a) ’) 
1 egend ( ‘ v_a_s ’ , ’ v_b_s ’ , ’ v_c_s ’ , ’ Locati on ’ , ’ Southwest ’ ) 
ylabel  (‘Voltage  (V)’);  xlabel  (‘Time  (s)’) 
xl i m ( [max (ti mev) - . 1  max ( t i mev) ] ) 

subplot(2,l,2) 

plot(timev,ias,timev,ibs,timev,ics,  ’Linewidth’  ,2)  ;  title(‘  (b)  ’) 
1 egend (‘i_a_s’ , ’i_b_s’ , ’i_c_s’ , ’Location’ , ’Southwest’) 
ylabel  (‘Current  (A)’);  xlabel  (‘Time  (s)’) 
xlim([max(timev)-.l  max(timev)]) 
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7,  qdO  Variables 


fig  =  fig+1; 

figure  (fig) 
subplot(3,l,l) 
plot(timev,i_qdos(:  ,1)) 
xlim([4.9  5]) 
legend(‘i_q’) 
subplot(3,l,2) 
plot(timev,i_qdosC:  ,2)) 
xlim([4.9  5]) 
legend(‘i_d’) 
subplot(3,l,3) 
plot(timev,i_qdos(:  ,3)) 
xlim([4.9  5]) 
legend(‘i_0’) 

iqs=i_qdos(: ,1) ; 

ids=i_qdos(: ,2) ; 

iOs=i_qdos( : , 3) ; 

theta_r=i_qdos ( : , 4) ; 

i qd_angl e=angl e (i qs-li *i ds) *180/pi ; 
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8,  Velocity  and  Position  (Electrical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,wre_data, ’Linewidth’ ,2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Angular  Velocity  (rad/s)’) 
xlim([4.9  5]);  %ylim([-125  125]) 

%1 egend ( ‘ \omega_r_e ’ ) 

subplot(2,l,2) 

plot(timev, theta_elec, ’ r ’ , ’ Li newi dth ’ , 2) 
title(‘(b)’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Electrical  Position  (rad)’) 
xlim([4.9  5]) 

%1 egend ( ‘ \theta_r_e ’ ) 

%ylim([0,2*pi]) 
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9,  Velocity  and  Position  (Mechanical) 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 
plot(timev,  V,  ’  Li  newi  dth  ’  ,  2) 
title(‘  (a)  ’) 

xlabel  (‘Time  (s)’);  ylabel  (‘Velocity  (m/s)’) 
xlim([4.9  5]); 

%1 egend ( ‘ Vel oci ty ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

pi  ot (ti  mev , z , ’ r ’ , ’ Li newi dth  ’  ,2) 
title(‘(b)’) 

xlim([4.9  5]);  ylim([0  2*pi/10]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Position  (m)’) 

%1  egend (‘Position’) 
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(a) 


10,  DC  Variables 


fig  =  fig+1; 
figure  (fig) 

pi ot (ti mev , vdc , ’Li newi dth ’ ,2) 

xlim([4.9  5]);  %ylim([0  11]) 

xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fig  =  fig+1; 

figure  (fig) 

plot(timev,  vdc,  ’  Li  newi  dth  ’ ,  2)  ;  %title(  ‘  (a)  ’) 

%xlim([0.5  0.6]);  %%ylim([0  11]) 
xlabel  (‘Time  (s)’);  ylabel  (‘Voltage  (V)’) 

fi  g=fi  g+1 
figure  (fig) 

plot(timev,idcabc(:  ,1)  ,timev,idcabc(:  ,2)  ,timev,idcabc(:  ,  3))  ; 
%title(‘(b)’) 

legend(‘i_a  dc’,’i_b  dc’,’i_c  dc’ , ’Location’ , ’Best’) 
xlim([2.34  2.367]);  %%ylim([-l  17]) 
xlabel  (‘Time(s) ’)  ;  ylabel  (‘Current  (A)’) 
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%  %%  Lambda  Harmonics 
%  fig  =  fig+1; 

% 

%  figure(fig) 

%  plot(funtime,0.1549*(fund+harm3+harm5+harm7) , ’Linewidth’ ,2) 
%  xlim([0  max(funtime)] ) 

%  %legend(‘\lambda_P_M’) 

%  xlabel(‘Time  (arbitrary)');  ylabel  (‘Magnetic  Flux  (wb)’) 

fig  = 

10 
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Current  (A)  Voltage  (V) 


Time  (s) 
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11,  Input  Torque  and  Force 


fig  =  fig+1; 

figure  (fig) 
subplot(2,l,l) 

plot(timev,input_torque, ’ Linewidth’ ,2) 

titleCCa)’) 

xlim([4.9  5]); 

%legend(‘\tau_i_n’ , ’Location’ , ’Southeast’) 
ylabel (‘Torque  (N\cdotm) ’ ) 
xlabel  (‘Time  (s)’) 

subplot(2,l,2) 

plot(timev, i nput Force, ’ r ’ , ’ Li newi dth ’ , 2) 

title(‘(b)’) 

xlim([4.9  5]) 

%legend(‘F_i_n’ , ’Location’ , ’Southeast’) 
ylabel  (‘ Force  (N)’) 
xlabel  (‘Time  (s)’) 
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